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Abstract 

A better grasp of the physical foundations of life is necessary before we can un- 
derstand the processes occurring inside a living cell. In his physical theory of the 
cell, American physiologist Gilbert Ling introduced an important notion of the rest- 
ing state of the cell. He describes this state as an independent stable thermodynamic 
state of a living substance in which it has stored all the energy it needs to perform all 
kinds of biological work. This state is characterised by lower entropy of the system 
than in an active state. The main contribution to this reduction in entropy is made 
by the cellular water (the dominant component with a concentration of 14 M) which 
remains in a bound quasi-crystallised state in a resting cell. When the cell becomes 
active the water gets desorbed and the system's entropy goes up sharply while the free 
energy of the system decreases as it is used up for biological work. However, Ling's 
approach is primarily qualitative in terms of thermodynamics and it needs to be char- 
acterised more specifically. To this end, we propose a new thermodynamic approach 
to studying Ling's model of the living cell (Ling's cell), the centrepiece of which is the 
non-ergodicity property which has recently been proved for a wide range of systems in 
statistical mechanics [7]. In many ways this new thermodynamics overlaps with the 
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standard quasi-stationary thermodynamics and is therefore compatible with the prin- 
ciples of the Ling cell, however a number of new specific results take into account the 
existence of several non-trivial motion integrals communicating with each other, whose 
existence follows from the non-ergodicity of the system (Ling's cell). These results 
allowed us to develop general thermodynamic approaches to explaining some of the 
well-known physiological phenomena, which can be used for further physical analysis 
of these phenomena using specific physical models. 
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1 Introduction 



The living state of a substance has always attracted the attention of physicists. And indeed, 
only a clear understanding of the thermodynamic characteristics of the living matter can 
give us insight into the processes that occur in living organisms. However, despite the fact 
that there is a lot of interest in this problem, this field can't be said to be developing by 
leaps and bounds. 

Before proceeding with our analysis we have to establish its boundaries and conditions. 
There are two very different approaches to the thermodynamics of living systems: one based 
on thermodynamics of equilibrium process [1,2,3,4] and the other on thermodynamics of 
non-equilibrium processes [5]. Shroedinger's research [6] differs from both of these. 

An obscure Russian scientist of Hungarian descent Ervin Bauer [4] was probably the first 
to suggest that the living state should be regarded as an unstable equilibrium. Treating the 
physical state of living substance in this way allowed him to draw a number of interesting 
conclusions and generalizations, but on the whole the most part of research was purely 
theoretical. 

According to Ling [1,2,3], whose position is of special interest to us, the minimal cell in the 
physical sense is a complex comprising protein in an unfolded configuration and water with 
ions , water and ions. The most significant characteristic of this complex is the state of water 
in it; it is absorbed by the protein in the form of a multi-layered structure that surrounds 
it along the entire length of the polypeptide. This 'coat' consisting of water molecules is 
stabilized by hydrogen links that are stronger than the hydrogen links in volumetric water. 
This increase in strength is a result of an increase in the dipole moment of the water molecules 
under the influence of other dipoles that are stronger than water such as the functional groups 
in the peptide link bound (NH and CO). The polarization of water molecules explains both 
their strong binding by the polypeptide frame of the protein and the multi-layered absorption 
of water on the surface of the unfolded protein. According to Ling, almost all the water in 
a the cell is in a bound state. Because, in terms of the number of molecules, water is the 
most abundant compound inside the cell, its transition into a quasi-crystal state results in 
a significant fall in the entropy of the cell. This lower entropy is what brings about the rise 
in the amount of free energy of the resting living cell. 

According to Ling [1,2,3], whose position is of special interest to us, the minimal cell 
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in the physical sense is a complex comprising protein in an unfolded configuration, water 
and ions. The most significant characteristic of this complex is the state of water in it; it is 
absorbed by the protein in the form of a multi-layered structure that surrounds it along the 
entire length of the polypeptide. This 'coat' consisting of water molecules is stabilized by 
hydrogen links that are stronger than the hydrogen links in volumetric water. This increase 
in strength is a result of an increase in the dipole moment of the water molecules under the 
influence of other dipoles that are stronger than water such as the functional groups in the 
peptide bound (NH and CO). The polarization of water molecules explains both their strong 
binding by the polypeptide frame of the protein and the multi-layered absorption of water 
on the surface of the unfolded protein. According to Ling, almost all water in the cell is in 
a bound state. Because, in terms of the number of molecules, water is the most abundant 
compound inside the cell, its transition into a quasi-crystal state results in a significant fall 
in the entropy of the cell. This lower entropy is what brings about the rise in the amount of 
free energy of the resting living cell. 

The introduction of the concept of a resting state is one of Ling's achievements. It's 
this state that is used as the reference point for all the physical and chemical processes 
that take place inside the cell. When a cell is activated by an external stimulant or some 
other signal, it changes its state from resting to active. The active state is characterized by 
the disintegration of the water-protein-ions complex. The bound water breaks free and the 
system's entropy increases. The free energy of the resting state is released and is used up 
for all kinds of biological work. This is the Ling model of the living cell (Ling's cell) that 
will be the focus of our analysis. 

According to Ling, a cell can remain resting without exchanging energy or substances with 
the external environment. The cell just maintains diffusion equilibrium with the. This view 
directly contradicts Prigozhin's approach according to which a living cell can only maintain 
its organization as long as it keeps exchanging substance and energy with the environment 
on a continuous basis. In other words, a cell can be likened to a burning candle flame; the 
flame will remain 'alive' only as long as there is sufficient supply of fuel and oxidizer. 

Ling believes that such understanding of the thermodynamics of life is completely inade- 
quate in the case of a living cell. His calculations demonstrate that if a continuous inflow of 
energy was really necessary for the experimentally observed exchange of Na+ ions between 
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a resting cell and the environment (as is postulated in the traditional mechanism), the cell 
would simply be incapable to produce the necessary amounts of energy [3] and therefore the 
universally adopted model of ion transport contradicts the energy preservation law. Ling's 
other argument proceeds as follows; if a living cell and a burning candle were to be frozen to 
the temperature of liquid nitrogen, both the flame and the life in the cell will 'go out', but if 
they're heated back to room temperature, the flame won't start burning again, but the life 
processes in the cell will resume. 

The contradictions between the thermodynamic approaches to the phenomenon of life 
are so pronounced that the need for further research in this field is self-evident. The purpose 
of this paper is the demonstrate that the property of non-ergodicity that has recently been 
proved for a large number of systems in statistical mechanics can help better understand the 
nature of the resting state of a Ling's cell and supports his understanding of the living cell's 
thermodynamics. 

One feature of this approach is that it suggests that the resting state of a Ling's cell 
should be considered to be a non-equilibrium stationary state, whose existence is a direct 
consequence of the non-ergodicity property that we postulate for Ling's cells. In this ap- 
proach thermodynamics of non-equilibrium stationary states must be constructed (analogous 
to the standard quasi-stationary thermodynamics) to explain why biological work becomes 
possible in the context of the proposed approach (biological work here means any changes 
in the cell that have a biological significance and that use up energy, for instance muscle 
contraction). There is no real contradiction between Ling's stationary resting state and the 
obviously continuous metabolism necessary to maintain life, because a real cell constantly 
changes its state from active to resting and back. Metabolism and energy are needed to go 
back to a resting state rather than maintain it. 

Non-ergodicity means that there exist non-trivial first integrals of the system (i.e. there 
are values that are invariant under the Heisenberg and Hamilton motion equations). Be- 
cause by definition we consider these first integrals to be real, then, in the case of quantum 
mechanics, they must be represented by self-adjoint operators and be experimentally observ- 
able values. It's then only natural to ask how come they can only be observed in biological 
systems (as is shown below) as well as in liquid helium and superconductors, systems that 
are as far from biological as can be? When answering this question we come across a certain 
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mathematical similarity between the super-fluid state of helium and the resting state of the 
Ling cell, and this similarity helps us better understand the physics of the living state. 

The main result of this paper is that by looking at the Ling cell as a non-ergodic system 
we were able to propose a common physical mechanism for various physiological phenomena, 
which were previously explained with the help of separate mechanisms barely related to each 
other. The main goal of physics in physiology must be to find out the thermodynamic nature 
of an active living cell, i.e. the source of all the manifestations of life. With this goal in 
mind we only discuss some of the characteristics of the Ling cell and provide only a most 
general physical description for them. We demonstrate, for example, that when a Ling's cell 
is activated it emits heat rather than absorbs it. When we consider the properties of the 
physical model we use, it becomes perfectly clear why it is that unfolded proteins that make 
up the structural foundation of a resting Ling's cell, begin to fold when it goes active, why 
potassium ions exit the cell into the environment and why an active cell changes its size 
(usually it shrinks) and what makes a cell dead. It is not our goal in this paper to compare 
the results we obtain for the Ling's cell with the properties of a real living cell. 

2 Non-ergodicity of statistical mechanics systems. 

In this section we first give the definition of ergodicity in statistical mechanics. We then 
formulate the main result of [7] and its classical analogue. After that we demonstrate how 
non-ergodicity follows from this result on the classical level. 

Definition. Suppose that a quantum system is described by Hamiltonian H and some set 
of self-adjoint integrals of motion K 1: ...,Ki commuting to each other. We say that the system 
is ergodic with respect to the integrals Ki,...,Ki, if each dynamical variable commuting with 
H, Ki,...,Ki, is their function, or in other words the joint spectrum of operators H, Ki,..., 
K\ is simple. 

The classical version of this definition can be easily obtained by replacing the word 
commutator with a Poisson bracket. 

In our case of Bose gas, apart from the Hamiltonian, there are the following trivial 
commuting first integrals: impulse P and the number of particles N. We will be looking at 
ergodicity relative to this set of first integrals. 
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Let us now formulate the main result of our work [7]. Let ty(x) and be secondary 

quantized wave functions of Bose gas that satisfy the canonical commutative relations 



[W(a;), = [y + (x),^ + (x')} = 0, 

Mx),* + (x')]=5(x-x>), (1) 

where the brackets designate a commutator and 5(x — x')- means Dirak 5 — function of a 
vector argument. An algebra generated by secondary-quantized wave functions is called the 
algebra of canonical commutative relations. The main result of [7] can be formulated as the 
following theorem. 

Theorem 1. There exist a linear functional (•) on the algebra of canonical commutative 
relations (within a formal perturbation theory applied to an interaction constant) such that: 

a) It is stationary, i.e. ([H, v]) = 0, where H is a Hamiltonian of the system and v is an 
arbitrary element of the canonical commutative relations algebra. 

b) (-) is translation-invariant. 

c) )(-) Commutates with the operator of the number of particles N, i.e. ([N, v]) = 0, 
where N is the operator of the number of particles in the system and v is an arbitrary 
element of the algebra of commutative relations. 

d) Satisfy to the week cluster property, i.e. 

lim / (fy ± (t,x l + 5 1 e l a)...^ ± (t,x n + 5 n e 1 a)}f(x 1 ,...,x n )d 3 x 1 ...d 3 x n 

\a\^oo J 

= J (^t,x il )...* ± (t,x ik ))(V ± (t : x ik )...y ± (t : x in )) 

xf(x 1 ,...,x n )d 3 x 1 ...d 3 x n , (2) 

where Si G {1,0}, i = l,2...n and 

i\ < i 2 < ... < ik, 

ik+l < ik+2 < ■■■ < in, 

{h,i 2 , ...,ifc} = {i = l,2...n\5i = 0} ^ 0, 
4+2, in} = {i = l, 2...n\5i = 1} ^ 0. (3) 

f(xi, ...,£„) - is a test function (i.e. it is sufficiently smooth and it decreases sufficiently 
fast at infinity with all its derivatives), t\ is a unit vector parallel to the x axis. 



e)(-) is not describe the Gibbs distribution. 

Before proving the non-ergodicity of the system from this result, we will formulate its 
classical analogue. In classical statistical mechanics a system containing N particles N — > 
+00, is described by a distribution function of the form p((xi,pi), (xn,Pn)), where (xi,Pi) 
are the coordinates and impulse of the z-th particle % = 1, N. The distribution function p is 
assumed to be symmetrical relative to the permutation of its arguments (xi,pi), (xn,Pn)- 

Also partial distribution functions are introduced: 

Po = 1, 

pi{x\) = V J dx 2 dp2...dx N dp N p(( y x 1 ,p 1 )...(x N ,p N )), 

P2((x 1 ,pi),(x2,p 2 )) = V 2 J dx 3 dp 3 ...dx N dp N p((xi,p 1 )...(x Nl p N )) 1 
(4) 

V is a volume of the system. Partial distribution functions p = 1, pi, p n ... satisfy the 
following consistency condition: 

Pni&uPi), (x n ,p n )) = 

= V m / dx m+ idp m+ i...dx n+m dp n+m p n ((xi,pi), (x n+m ,p n+m )). (5) 



If the number of particles in the system is variable, then in a similar way to (j3J), we can also 
define partial distribution fucntions. They will also satisfy the consistency conditions if N 
approaches infinity. If we take the main result of [7] and pass it through the limit we get the 
following theorem. 

Theorem 1'. There exists such a state p (not necessarily positively defined and possibly 
with a varied number of particles) of a gas consisting of particles with weak interaction, that 

a) This state is stationary, i.e. it doesn't change with time due to its motion equations. 

b) This state is translation- invariant. 

c) This state satisfies the weak cluster property, i.e. Wn = 1, 00 and Wk = 1, n 



lim p n ((xi,p{), (x k ,p k ), (x k +i + le^pk+x), (x n + le x ,p n )) = 

\l\— s>oo 

= p k ((xi,pi), (x k ,p k ))p n -k((x k+1 ,p k+1 ), (X n ,p n )) (6) 

d) The state p cannot be described with a Gibbs distribution. 

8 



Let us finally demonstrate how this result leads to the non-ergodicity of our system. 
To simply things we'll only be considering the case of classical mechanics. For a detailed 
discussion of the quantum version of this result see [7] . 

Let us suppose that our system is ergodic, i.e. that its energy and impulse form a maximal 
set of first integrals in involution. Two dynamic variables are said to be in involution if 
their Poisson bracket equals zero. From here on for the sake of brevity, we'll be simply 
talking about energy rather than energy, impulse and the number of particles. Let / be a 
distribution function for our system in a phase space that corresponds to the state we talked 
about before. Then because of the system's assumed ergodicity it must be a function of only 
energy / = f(E r ) (F - is a point in the phase space) and can be presented as a superposition 
of micro-canonical distributions: 



The summation in this context is understood in the broadest possible sense and can be 
infinite (in other words it can be replaced with integration). 

Let 1 be a sufficiently large but finite subsystem of our system and let 2 - be a subsystem 
obtained from 1 by shifting it along the x axis by a distance of I. Let 12 - be a union of 
systems 1 and 2. Let Ti, T 2 and Ti 2 designate points in the phase space for our systems 1, 
2 and 12 respectively. Let /i(Fi), 72^2) and /^(r^) be distribution functions for systems 
1, 2 and 12 respectively. Then using the same method that is used for obtaining a canonical 
distribution from a microcanonical one, for |/| = 00 we can find that: 



Here d a > - is a weight multiplier, Er 1 and Er 2 - are the energies of subsystems 1 and 
2, T a - is the temperature corresponding to energy E a , and Z a - is the statistical sum: 




(7) 




(8) 




(9) 



But from the weak cluster property it follows that: 



fl2 — flf2- 



(10) 



9 



Therefore it is follows from (jSJ) that all c a = except one and: 

f(E)=c8(E-E ), (11) 

(for some c and E Q ), i.e. the whole system is described by a microcanonical distribution 
while any of its sufficiently large but finite subsystems is described by a canonical Gibbs 
distribution, which contradicts the fact that our state is not a Gibbs state. This contradiction 
proves that our system is non-ergodic. 

3 Discussion of Non-ergodicity 

Because not enough research has bee conducted into the specific features of the dynamics 
of statistical mechanics systems connected with non-ergodic property, the discussion in this 
section will be informal and hypothetical in places. 

Boltzmann's hypothesis. Boltzmann hypothesis (1871) postulates that in statistical 
mechanical systems the dynamic variables that are average in time are equal to the dynamic 
variables that are average across the (micro) canonical ensemble. We will only consider 
the quantum mechanics case. The discussion presented here was taken from [7]. Before we 
can formulate the Boltzmann hypothesis for a Bose gas with weak pair interaction, some 
designations need to be introduced. Let (a)^^ be the Gibbs state corresponding to the 
inverse temperature (3, the system velocity v and its chemical potential fx, i.e. 

<«W = ^tr(ae-^-^)), (12) 

where a is an arbitrary element from the canonical commutative relations algebra, H is the 
Hamiltonian of the system, iV is the operator of the number of particles in the system, P is 
the operator of the system impulse and Z is the so called large statistical sum, i.e. 

Z^ = tr(e-^ N+ ^). (13) 

Let Vq be a space of all linear functionals on the algebra of canonical commutative 
relations, each of which is a superposition of Gibbs states. 

Let (•) be is a translation-invariant functional on the algebra of canonical commutative 
relations such that Vt G R (e itH \-)e~ itH ) is a well defined linear functional. The Boltzmann 
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hypothesis states that in this case there will be a linear functional (•)' G V G such that for 
any element a from the algebra of commutative relations: 

T 



Let us demonstrate that if formulated in this way the Boltzmann hypothesis is wrong. 
Theorem 1 states that there exists a translation invariant stationary linear functional on the 
algebra of canonical commutative relations and that said functional is not a superposition of 
Gibbs functions. However, this functional is built only in the form of a formal power series 
on the interaction constant and nothing is known about whether this series has a finite sum 
or not. [7], however, demonstrates that there is a way around the issue of whether this series 
converges or not and that there exists a stationary translation invariant functional on the 
algebra of canonical commutative relations (-) s that does not belong to Vq. If this composite 
function is substituted for (•) in the left side of equation (|T4j) . then the left side itself will 
equal (-) s and from equation f|T4|) it will then follow that (-) s is superposition of Gibbs states. 
This contradiction shows that the Boltzmann hypothesis does not hold. 

Discussion of the relation between a system's behaviour on its boundaries 
and its tendency towards thermodynamic equilibrium. Non-ergodicity means that 
Bose gas (with weak pair interaction) will not tend towards thermodynamic equilibrium in 
an infinite volume. Consequently to prove the tendency towards thermodynamic equilibrium 
we have to take into account its behaviour on the boundaries. In fact if the role played by 
the boundaries is disregarded we may end up with an infinite system. 

In [7] the role of boundaries is illustrated with Bogolubov's derivation of Boltzmann's ki- 
netic equation. Bogolubov's programme for deriving equations considers an infinite strength- 
ening chain of equations for partial distribution functions (BBGKI chain) equivalent to 
Luiville's equations, and then closed kinetic equations are obtained by selecting the suitable 
conditions for breaking the chain (correlation breaking conditions ). The first member of 
this chain has the form: 



Here t is time, m is the mass of the particles, and V is the potential of pair interaction. 




(14) 



o 



d_ 




(15) 



n 



An integral over the coordinates in the last term in the equation on the left ( I15p may be 
transformed into an integral along the boundary of a three-dimensional space using the 
Gauss theorem, and usually this fact is used as the basis for assuming that this integral 
equals to zero. However in [7] it is demonstrated if the last term on the left side of ( TT5l) 
is not discarded and Bogolubov's strategy for deriving a Boltzmann equation is followed 
(if translation- invariant distributions are considered from the start), then this term will 
eventually cancel out with the collision integral and the system will have no kinetic evolution. 
Even though, it is possible that this effect is purely mathematical, it does make one wonder 
about the role that the system's behaviour on its boundaries plays in proving its tendency 
towards thermodynamic equilibrium. 

Although it is now clear that the system's behaviour on the boundary must be taken 
into account when proving its tendency towards thermodynamic equilibrium, the following 
question now arises: if the linear sizes of the system equal ~ L and tend to infinity, then the 
volume of the system will b ehave like L 3 while the surface area of its boundary will equal 
L 2 which makes one think that the influence of the boundaries must be very insignificant. 

A similar situation can be observed in the theory of second order phase transitions, for 
instance in the theory of ferromagnetism. If, for instance, you take a two dimensional Izing 
model at a temperature below the critical level, and assume that all the spins are directed 
upwards along the boundaries of the crystal, then Peierls estimates (see [8]) show that in the 
thermodynamic limit there will be a non-zero magnetization. Second order phase transitions 
are often related to breaking of some symmetry. In the case of the Izing model, the symmetry 
that is broken is the symmetry relative to the change in the direction of all the spins in the 
grid. In the case of kinetic evolution there is also as symmetry that gets broken, it's the 
symmetry relative to the time sign (Boltzmann's H-theorem). We will have a more detailed 
discussion of this similarity later on. 

The necessity to take into account a system's behaviour on its boundaries when proving 
the system's tendency towards thermodynamic equilibrium raises an interesting philosophical 
question. If a system's behaviour on its boundaries must be taken into account when proving 
the system's tendency towards thermodynamic equilibrium, then it means we have to know 
about the behaviour of the particles outside the system's boundaries because those influence 
the behaviour of the particles immediately inside the system's boundaries, in other words 
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we have to know how the environment that the system is a part of behaves. This in turn 
means that we have to know the behaviour of the larger environment around the immediate 
environment of the system and so on and so forth. So how do we escape from this infinite 
regress? 

To answer this question we will consider the following example. It is a known fact 
that there are black holes in our universe. Black holes are known to emit energy but for 
this phenomenon to be explained certain boundary conditions have to be assumed, which 
essentially are based on the fact that a black hole is formed as a result of the collapse of 
a mass body. According to one theorem of Penrose, this means that at one time in the 
past there existed a certain feature of the space time. The surface of a black hole, the 
so called event horizon, is the boundary at which the dependence of the particles inside 
the system on particles outside the system is broken. In fact due to the structure of light 
cones, information on the event horizon can only fall into the black hole, while the Hawking 
radiation is in equilibrium and therefore carries only the bare minimum of information. 

On the other hand one of Penrose-Hawking's theorems [9] (whose conditions are con- 
firmed by reliable astronomical observations) postulates that at some time in the past there 
was a singularity of space-time which can be equated to the Big Bang. At the same time 
Boltzmann also posited (on the basis of his i^-theorem) that there had to have been a mo- 
ment in time which can be equated to the beginning of the universe. In the theorems of 
Penrose and Hawking there appears entropy (equal to the surface area of the event horizon) , 
which increase and is related to the temperature of the Hawking radiation through a known 
thermodynamic relation. It follows from all of this that there is a much deeper link between 
the 'proof of Boltzmann's Big Bang and Penrose — Hawking's singularity than a simple 
mathematical similarity. 

It is also worth noting that due to the arguments presented above, the non-ergodicity 
of Bose gas with weak pair interaction can be interpreted as the impossibility of deriving 
macroscopic dynamics from microscopic dynamics. 

The non-ergodicity theorem, as follows from its proof (see [7]), is applicable to a broad 
range of systems in statistical mechanics and we make the assumption that it can also be 
applied to the Ling's cell. Therefore, the general results obtained for all non-ergodic systems 
can also be applied in this case. 
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A more detailed study of how the properties of a system's surface influence the process 
whereby it achieves thermodynamic equilibrium can be helpful in improving our understand- 
ing of living cells and the role of the cell's surface in the physiological processes inside it. 

Analogy with second order phase transitions. Let's discuss in more detailed the 
aforementioned analogy between the broking of the symmetry relative to the time sign and 
second order phase transitions. We'll do this by discussing the Green-Cubo formula that 
describes the linear response of a system to small perturbations. Let there be n fluctuating 
variables in the system xi,...,x n and let Xi,...,x n be corresponding quantum mechanics 
operators. Let us assume that the fluctuations are sufficiently small and that the system's 
entropy can be expressed with the formula 

S = Sq — PijXiXj, (16) 

where is some positively defined matrix (here and through the end of this section indexes 
repeated twice mean summation). Let X±, ...,X n be thermodynamically conjugate values to 
X\ , . . . , x n i.e. 

= ( 17 ) 

ox k 

The equation describing the relaxation of this system look like this: 

ii = -jikXh, (18) 

where 7^ - are the so called kinetic coefficients. This coefficients satisfy Onsanger's 
symmetry principle: 

lik = Iki- (19) 

If the system is affected by external forces f k k — 1, n, in other words if the system's 
Hamiltonian looks like this 

H = H + f k x k , (20) 

(H is the Hamiltonian of an unperturbed system), then we can determine the generalized 
sensitivity of the system 0^(7-) in such a way that 

+00 

x k (t) = J a ki (t - r)/i(r)dr. (21) 

—00 
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Let's use a ik (u) to designate the Fourier transform of the generalized sensitivity 

atikiu) = J e iut a ik {t)dt. (22) 
The generalized sensitivity can be found from the Green — Cubo formula: 

a ik (t) = -i0(t)<[£i(t),£ fc (O)]), (23) 

where 6(t) is a Heaviside step function that equals 1 when its argument is positive and zero 
if its negative and (■) is the averaging over the equilibrium state of an unperturbed system. 
It has to be noted, however, that when the Green — Cubo formula is applied the system is 
assumed to be in equilibrium with t = —00. If we assume, on the contrary, that the system 
is in equilibrium with t — 00, then by transforming the time in the Green Cubo — formula, 
we will get the following formula for the generalised sensitivity 

a' ik (t) = te(-t)([x l (t),x k (0)]). (24) 

It can be easily seen that 

a' ik {u) = (a ik {cu))\ (25) 

where the asterisk means complex conjugation. If, however, we start with an approximation 
in which relaxation equations ffl8|) hold, then it can be shown that there exists the following 
relations between the matrices a, (3 and 7: 

a ik {u) = ^{f3 ik - iw7^ 1 ) -1 . (26) 

It can be seen from this formula that if t is replaced with —t and vice versa in the Green - 
Cubo formula, then j ik will become —7^ and will no longer describe kinetic evolution. 

So, how can we choose between the requirement for equilibrium with t = —00 and 
with t = +00, in other words how can we determine whether the equation with a retarded 
commutator ( 1231) or the equation with a advanced commutator ( 1261) should be used for the 
generalized sensitivity of the system? 

It should be noted that a similar problem arises in classical electrodynamics when working 
out the formula for retarded potentials. A vector potential A 11 , fi = 0, 1, ...,3 in classical 
electrodynamics satisfies the heterogeneous wave equation: 

DA* = J M , (27) 
15 



where 

□ = ? ? ^_ f28) 

(ax ) 2 (5s 1 ) 2 (ax 2 ) 2 (ax 2 ) 3 ' 1 ; 

and J M is a four- vector of current, satisfying the continuity equation: 

dx ^ ~^ dx 1 ^ ~^ dx 2 ^ ~^ dx 3 '^ ^ ^ ^ 

Where Xq is time, and x 1; X2,X3 are the spacial Euclidean coordinates. 
The solution to equation (1251) usually has the form: 

A»(x) = J D(x -y)J»(x)d 4 x, (30) 

where D(x) is the fundamental solution to the wave equation: 

DD(x) = 5(x), (31) 

with a support in the upper light cone: 

suppD(x) C V + = {x|x > 0, x 2 > 0}. (32) 

There exists a single fundamental solution that satisfies this property. Formula (1301) is called 
the formula of retarded potentials. The similarity to the Green — Cubo formula becomes 
even closer if we remember that the vacuum mean of the commutators of the secondary 
quantum vector potential in electrodynamics is presented in the following form 

(0\{A»(x),A»(y)]\0)6(x -y )= ig ^D(x-y), (33) 

where g^ v is a metric tensor: 

gT = 0, if fj, ± V, 

g 00 = l, g« = -l, 2 = 1,2,3. (34) 
If in (12*51) we switch to the Fourier transform in time, then we will get: 

{ul + Wf + Wf- + Wf }A ^^ = - J "'"- f) - (35) 

The solution to this equation, corresponding to the formula of delayed potentials, will have 
the form: 

A»(u,x) = -^J ^—^r(u;,y)d 3 y. (36) 
16 



If Sr is a sphere with a radius R and its center at zero and R — > oo, then A^(u,x) on this 
sphere will satisfy the Sommerfeld radiation condition: 

-«M>,aO = 0(-^), (37) 

where J| means differentiation along the external normal to sphere Sr. By applying a inverse 
Fourier transform to this formula, we will find the condition that must be met by Sr, to 
ensure that a retarded commutator is selected in the formula for retarded potentials: 

^ W + ^, (x) = 0( i_). (38) 
The selection of a advanced commutator would result in the following formula 

Thus we can see that a advanced or a retarded commutator in Green — Cubo formulas can 
be selected (which is necessary to specify the direction of time), at least in electrodynamics, 
by applying appropriate boundary conditions on the infinitely removed boundary of a three 
dimensional space. This approach is in full accord with our ideas about the role that the 
boundaries of a system play in proving its tendency towards thermodynamic equilibrium. It 
should also be noted that boundary condition (1381) imposes a limit on the flow of information 
(this interpretation is also possible): information can only go out of the system (flow outside). 
This restriction on the direction of the flow of information is realized, for instance, on the 
event horizon of a black hole due to the structure of the light cones on it, as has already 
been noted above. 

Let us go back to our analogy with ferromagnetism (Izing model). We said earlier that 
non-zero magnetization can be obtained by changing the direction of the spins on the bound- 
ary of the grid. However, non-zero magnetization can also be achieved in another way. We 
can introduce a fictitious infinitesimally small external field, which is expressed by adding 
the following expression to the system's Hamiltonian 

hJ2<r p , (40) 
p 

where the summation is done for all the nodes p of the grid, and a p is the spin in the node 
of the grid with a p G { — 1,1}. Then the system passes to an unlimited volume and after 
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that h starts approaching zero, however it always remains positive or negative. The means 
calculated in this way are called quasi-means (Bogolubov's quasi- means) . It turns out that 
non-zero magnetization occurs in this case as well (and it's equal to the magnetization that 
results when the spins on the boundaries are given a certain direction). 

It begs the question whether it is possible to make the right choice between a advanced 
and retarded commutator method when working out a Green — Cubo formula, by analogy 
with Bogolubov's quasi-means? The answer to this question is positive and here is why. Let 
H(t) be a Hamiltonian of the system dependent on time: 

H(t) = H + Vf(t), (41) 

Let f(t) be an external and infinitesimally small classical force changing with time. The von 
Neumann equation for the density matrix will look like this : 

j f p(t) = -i[H(t),p(t)} = -i[H ,p(t)]-if(t)[V,p(t)}. (42) 

If f(t) = 0, p = p is the density matrix of the equilibrium of an unperturbed system. 
Consequently if we're looking for p in the first order over /(£), then the following equation 
can be used for p: 

j t p(t) = -i[Ho,p(t)}-if(t)[V,p }. (43) 

The standard method for solving this equation when working out a Green — Cubo formula 
involves choosing a boundary condition p(t) — > p with t — > — oo and solving ( 14"3"|) by using 

the constant variation method. What we get then is: 

t 

p = p -i J e^ Ao ^[V,p ]e +t6o{t ~ T) f(r)dT. (44) 

— oo 

This expression for the density matrix immediately leads to a Green — Cubo formula with a 
retarded commutator. Instead of imposing initial conditions on the density matrix, we will 
introduce an infinitesimally small term ep(t), e > 0, e — > into the left side of (1431) . In other 
words instead of f|43|) we will be solving the following equation for the density matrix: 

lp(t)+ep + i[H ,p(t)} = -if(t)[V,p }. (45) 

Applying a Fourier transform to this equation we get: 

(-iw + e + i[H Q ))p(u) = -if(uj)[V,p l (46) 
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where [H ] designates the following super-operator: 

[H ]p := [H , p]. (47) 

Therefore it follows from (1461) that: 



PM = , - - , \V,p Q ]f{u). (48) 
-lu + e + i[H \ 

Since a Fourier transform transforms convolution into product we can find: 

+00 

p(t)= J a(t-r)f(T)dr, (49) 

—00 

where 



-00 



a(t) = i- / e**- t -- w -[V,p }du. (50) 

^ J -ico + e + i[H ] 

—00 

The super-operator [Hq] is the Liouvillian of an unperturbed system and is therefore self- 
adjoint (if we define the scalar product of the operators applicable in the space of states as 
(A, B) = tr(A*B)). Then, calculating the integral in (150]) through residue we will get: (150|) 



a(t) = -i6{t)e { - i[ilo] - £)t [V , po] = 6(t)e' £t e^ Ant [V, p )e i6ot (51) 

Substituting the expression we just found in ( 1491) for a(t) and assuming that e — > we will 
find: 

t 

P = P*-i I e- iA ^[V, Po \e^ {t -^f{r)dT, (52) 



-00 



Thus we've once again got expression ( )44l) for the density matrix. Thus we see that boundary 
conditions for a density matrix can be chosen using a method similar to Bogolubov's method 
of quasi means, i.e. by introducing an appropriate infinitesimally small term in the equa- 
tions describing the evolution of the density matrix through time. This should explain the 
similarity between spontaneous breaking of the symmetry in second order phase transitions 
and the breaking of the time symmetry in kinetic equations. 
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4 Thermodynamics of stationary non-equilibrium states 

In section 1 we said that according to the results obtained in [7], even the most realistic 
systems in statistical mechanics are non-ergodic. It should also be reminded at this stage, 
that a system is ergodic if its trivial first integrals (Hamiltonian, impulse, number of particles) 
form the maximum set of first integrals in involution. Hereinafter for the sake of brevity 
we will use the term energy to mean the energy and the impulse of a system. From the 
non-ergodicity of a system it follows that it must have stationary states that are not super- 
positions of Gibbs states. We will call the states of a system that are super-positions of 
Gibbs states equilibrium states. 

It therefore makes sense to build thermodynamics of non-equilibrium stationary states. 
This thermodynamics will really come in handy in the next section when we formulate the 
principle that the substance inside a living cell is in a non-equilibrium stationary state, which 
is the goal of this paper. This section describes how we build this thermodynamics. The 
equations we derive in this section will be applicable to any system in statistical mechanics, 
including the Ling's cell. For the sake of simplicity we will only consider the case of classical 
mechanics. 

Some basic concepts of classical (Hamiltonian) mechanics should be remembered here. 
The phase space T of a Hamiltonian system is a 2n-dimensional n — 1,2, ... smooth manifold 
in which for any two (good) functions f, g there is a function (/, g) defined on their Poisson 
brackets, which satisfies the following conditions: 

1) (fg, h) = f(g, h) + g(f, h) (Leibniz rule) 

2 ) (f,g) = ~(g,f) (Anti-symmetry) 

3) (/, (g, h) + (g, (h, /)) + (h, (f, g)) = (Jacobi identity) 

4) For any point x E T, and for any function /: grad/(x) ^ 0, there is a function g on T 
such that (f,g)(x) ^ (non-degeneracy). 

n is called the number of degrees of freedom of the system. 

An example of a phase space would be Y = {(pi, ...,p n ,qi, ■■■,q n ) \ Pi,q% € M} in which 
the Poisson bracket is defined as: 



According to the Darboux theorem (see [10]) local coordinates p±, ...,p n , qi, q n can always 




(53) 
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be introduces in which the Poisson bracket will look like (153 p . Such coordinates are called 
canonical or symplectic. 

The evolution of a Hamiltonian system is described by a one-parameter group of diffeo- 
morphisms (i.e. homeomorphisms which are smooth with their inverse), x h- > G T (x), r G R 
preserving the Poisson bracket (r is time). This one-parameter group of diffeomorphisms is 
called a phase flow. 

It can be shown that for any phase flow G T (locally) there is a function H on T, called a 
Hamiltonian, such that 



Where f T (x) := f(G T (x)). For more about Hamiltonian dynamics see ([10]). 

Thus, let's assume we have a Hamiltonian system in which in addition to Hamiltonian H 
there is also , k — 1, 2, 3... first integrals Ki, ...Kk in involution, i.e. all the matching Poisson 
brackets equal zero: Wi,j = 1, k (Ki, Kj) = 0. We want to describe a state of our system 
in which it is in equilibrium with its environment and in which its Hamiltonian H and the 
first integrals K±, Kk take on some specific values E, K[, K' k respectively. 

The requirement that the first integrals K\, Kk must be in involution can be explained 
as follows: we want to describe states in which the integrals of motion Ki, Kj, take on spe- 
cific values: E, K[, ■■■,K' k . But in quantum mechanics (see [11]) simultaneous measurability 
of the observables means that they are commutative. Now it becomes clear why K\, Kk 
must be in involution, if we take into account the fact that the commutator is the quantum 
mechanics counterpart of the Poisson bracket (see [12]). 

Let us first consider the way the above task is solved with k = 0. The state of a system 
in statistical mechanics is described by a distribution function p(x) on the phase space T, 
i.e. by a function that satisfies the following conditions 

a) Vx G Tp(x) > (Positive definiteness) , 

b) J p(x)dT x = 1 (Normalization). 

Here the phase volume is determined in the following way: in canonical coordinates it 
looks like dT = dpi...dp n dqi...dq n . It can be demonstrated that this definition is correct. The 
observed macroscopic value F of dynamic variable F(x) has the form: 



d 

d^ 



fr = (H,f T ) 



(54) 




(55) 
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The role of the distribution function p(x) is assigned to the so called micro-canonical Gibbs 
distribution: 



p(x) = c5(H(x) — E), (56) 

In which the constant c is chosen on the basis of the normalization condition. 

It's important to note here that we assumed that the system under consideration is in 
thermodynamic equilibrium with its environment. For this reason H(x) takes into account 
interaction with the particles in the environment and H(x) is a function of time and time is 
included in H(x) through the coordinates and impulses of the particles in the environment. 
This will be important further down the road. It has also be remembered, however, that in 
the thermodynamic limit these extra terms in H(x) that account for the interaction with 
the environment, can be disregarded. 

If the system's Hamiltonian is a function of time, H(t), then the system's evolution will be 
described by a co-cycle, i.e. Vti,t2 there are canonical transformations (i.e. transformation 
that preserves the Poisson bracket) x (-)■ G tlt t 2 (x) meets the co-cycle property: 

Any dynamic variable / is a function of time by definition ft(x) := f(Gt,o(x)). Now it has to 
be said that in actuality, the motion of the system is described by some trajectory through 
the phase space and for this reason the left side of (I55I) has to be interpreted as a mean over 
time, or 

T 

lim ^ J F(G Tfi (x))dr = cj F(x)S{H(x) - E)dT x . (58) 
o 

Equation fl58l) does not in any way mean that the system is ergodic, because from the 
very beginning we said that our system is a subsystem of a large system. 

A phase flow preserves its phase volume(see [10]). Therefore if is the system's distribution 
function, then p Consequently, over time the distribution function will evolve into: 

p t {x) = p(G 0)t (x)). (59) 

Function pt meats the conditions of the following differential equation: 

lp t = (p u H(t)). (60) 
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Now let us assume that the system's Hamiltonian H\ is a function of some parameter A. 
We will demonstrate that if A undergoes adiabatic changes the micro-canonical distribution 
remains the same. Let X(t) be some function of time and let e be some small positive number. 
Let \ £ (t) := X(et). Let a; be a point in the phase space such that H X ( £ t)(x) = E. Let [0, A] be 
a sufficiently long period of time such that a mean over this period is the same as the mean 
over the micro-canonical ensemble. Let us now select such a small e that will not change 
very much in this time period. Let's calculate 

H Xe(t) (G t , (x))\ t=A . (61) 

What we get is: 

H Mt) (G tfi (x))\ t=A = 

A 

= H ((x)) + J dt{f t d -^ + (H Xe{t) ,H Xs{t) )}(G t , (x)) = 
o 



= *„((*)) + A£S +0(e ») (62) 

Thus with a precision of up to e 2 canonical transformation G tt o transforms function H\ s ( t )(x) — 
E — A^- dH g^ {t) into H\ e (o-)(x) — E, i.e. transformation Gt,o transforms the micro-canonical 
distribution c5(H\^o)( x ) — E) into c5(H\ e ( t ){x) — E — A^ 9H ^ (t) ) with a precision of up to 
K The function \ £ (t) changes by a finite amount within an interval of order K Therefore, if 
we break down this time period into sub-intervals of duration A we will find that after time 
period ^ elapses the form of the micro-canonical distribution will change by a value of order 
e. If we have an e approaching zero, we can prove the above contention. 

Now we will consider the case when the system has k- first integrals in involution. By 
analogy with the micro-canonical distribution we will assume that in this case if the values 
of the parameters E, K[,...,K' k are fixed the system's distribution function will have the 
form: 

k 

p(x) = c8(H(x) -E)]J 5(K t (x) - Ki), (63) 

i=i 

wherec is a normalizing multiplier. First of all let us note that this definition is correct, i.e. 
the form of the distribution will not change if the integrals K ± , K k are replaced with the 
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functions K 



K k of these integrals and the Hamiltonian. In fact: 



k 



c8{H-E)Y[5{k i -k^ 



i=l 



k 




= Ci 



■8(H -E)Y[ 8(Ki - K-) 



(64) 



i=l 



where 



D(H,K u ...K k ) 

D(H,k 1 ,...k k ) 



(65) 



- is the Jacobian of the transformation of the variables H, Ki, ...K k into the variables 
H, K\,...Kk- This Jacobian is constant on the shared level surface of the functions H, Ki,...Kk, 
which proves the correctness of the definition. 

We will call distribution (1631) the generalized micro-canonical distribution. It has to be 
said that the generalized micro-canonical distribution can be derived from the requirement 
that the Von Neumann entropy must achieve its maximum if the energy value and the values 
of the first integrals Ki,... } K k , are fixed. In other words, the generalized micro-canonical 
distribution corresponds to a state of the system in which the information that is available 
about it (with the fixed values of the Hamiltonian and the first integrals Ki, is minimal. 

Using a method similar to that we used above, we can demonstrate that the micro- 
canonical distribution does not change its form as a result of quasi-static processes. The 
following formula has to be noted here: if the Hamiltonian H and the first integrals K\, K k 
are functions of parameter A, then in quasi stationary processes: 



where qs means that the derivative is calculated along the quasi-static process. We will now 
define the entropy corresponding to the micro-canonical distribution as: 




dK' 



(66) 




(67) 



Where W(E, K[, K' n ) is the so called statistical weight: 




i=l 
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Let's found out how correct this definition is, i.e. how will the entropy change if we move 
from the integrals H, Ki, to the integrals H, Ki, K^. What we get is: 



/n 
dTJ(H(x) - E)Y[5(ki(x) - K[ 
i=i 



(69) 



where c is the transformation Jacobian. The appearance of this multiplier next to the 
statistical weight will create an additional term in the entropy In c. It has to be noted here, 
that a similar problem also arises in standard thermodynamics (k = 0), where this parameter 
is simply discarded because it is assumed that it's influence in the thermodynamic limit is 
infinitesimally small. We will do the same here. It also has to be noted that instead of ( |63l) 
the following formula can be used for the micro-canonical distribution: 

k 

p{x) = cA(H(x) - E)Y[A(Ki(x) - K[), (70) 

i=i 

Where function A equals one in the small neighborhood of zero and equals zero outside 
this neighborhood. In the thermodynamic limit, this neighborhood will be so narrow that 
A-function can be replace with 5-function, and we will get back to ( 163]) . In light of the above 
notes, the statistical weight can be defined as: 

W(E,K[,...,K' n )= I dT x A(H(x) - E)Y[A(K t (x) - K^). (71) 
^ i=i 

It can be seen from this formula that the statistical weight can be interpreted as the number 
of microscopic states compatible with a given macroscopic state, or as the thermodynamic 
probability. 

Now let's assume our system is in equilibrium with the environment. Let's demonstrate 
that we can choose the first integrals Ki, K n in such a way that their values will remain 
the same as the system heats up. And indeed, the energy level E corresponds to the integral 
values K[(E), K' n (E). Let's replace the integrals K\, K n with the integrals K\, K n 
using the formula Ki = Ki — K'^H). It can be seen that these integrals remain constant as 
the system heats up. We then assume that these integrals also remain constant when the 
system temperature increases. What we then have is 



.3KL OK' 



v d\ >qs d\ 
From which we can derive: But because: 

,dKL ,dKL ,dKL .BE 



(72) 



8X Jqs v dX ,CJ ' ' OE OX 
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e + (^r^)\(-Kr)qs, (73) 



and (t^Oa = we have 



dK> 




dK[ 8Ki 



(74) 



dX 



dX Jqs dX ' 



It has to be noted that the entropy remains constant in a quasi-static process. Indeed, 
in a quasi-static process, the (generalized) micro-canonical distributions corresponding to 
the various parameters turn out to be bound, as was demonstrated above, by a canonical 
transformation and their entropy is the same because canonical transformations preserve the 
phase volume. One important result of this section is the conclusion that in the thermo- 
dynamics of stationary non-equilibrium states, the standard thermodynamic equations hold 
true. Let's derive the standard thermodynamic equation: 



where T is the temperature, P is the pressure, and V is the volume of the system. Put by 
definition that 

i = ^ (76) 

V dE K ' 

and demonstrate that T' = T, i.e. that it is the same as the ordinary temperature. Let 
E' be the energy of the environment around the system and S(E') its entropy. We will 
assume that our system is in equilibrium with its environment (but it's not in equilibrium 
with itself, being in a stationary non-equilibrium state.) The probability P(E, E') that the 
system's energy equals E, and the energy of the environment equals ~ e s(E)+s'(E') b ecause 
we earlier defined the system's entropy S(E) as the logarithm of thermodynamic probability. 
If the system is in equilibrium with its environment, then P(E, E') reaches its maximum 
(with an additional condition that E + E' = const), because equilibrium is the most probable 
state in this case. Consequently: 



dE = TdS - PdV, 



(75) 




(77) 



i.e. 



T' = T, 



(78) 



and that is exactly what we've been trying to prove. 
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Let us remind you that S = In W, where 

/n 
dTJ(E - H(x)) 11 5{K[{V) - Ki(x)). (79) 
i=i 

Here we directly specify the dependence of Ki on V . We just demonstrated that (f§)v = ^. 
Now all we've got left to do is to demonstrate that (fpr)js = ^-P- Hereinafter derivatives are 
taken with V = Vo for some Vo- What we get is: 

OS 
{ dV )E ~ 



i=l 



3=1 4=1, 4^j 



T7F / rfr :c {-^^^-if(x))n^-^)) + 

4=1 



Because of equation ( 1731) the second to last term on the right of the last equation ( 1801) equals 
zero. As far as the last term is concerned, the most reasonable assumption would be that 
^hrr\ ~ 1 or d Jt 3 }nX — 1 — const v ~ v ° , with V ~ Vn. Therefore the last term in (81) 

dKjiVo) dK'^Vo) Vo ' u v ' 

disappears in the thermodynamic limit. Thus the final result is: 



,dS_ J_ 
"dV' E ~ W I ' ' OY 



dT x ^^5'(E - H(x)) J] KK[ - K(x)) 



4=1 
n 



J ' j = l 



But — ( 8 gy ), as was already demonstrated earlier, is the speed at which the energy changes 
in an adiabatic process, in which the changing parameter is the volume, i.e. pressure. Thus 

But energy is an extensive value while pressure is an intensive one. For this reason in the 
thermodynamic limit |^ = 0. Finally, remember that = ^ we get: 

,dS x 1 . . 

( W h = f P, (83) 
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This is what we've been trying to prove. U 

sing the standard method and applying the Legendre transformation we can introduce 
thermodynamic functions. The free energy F is defined as: 

F := E — TS. (84) 

Using f!75|) we can find: 

dF = -SdT - PdV. (85) 
Similarly we can define the thermodynamic potential <3>: 

$:=F + PV. (86) 

What we get is: 

d$ = SdT + VdP. (87) 
If the number of particles in the system iV is variable then (1751) should be changed to: 

dE = TdS - PdV + fidN. (88) 
Using the theorem of small additions we can find that: 

dF = -SdT - PdV + fidN, 

d$ = -SdT + VdP + fidN. (89) 

Integrating the last equation with constant V and T between zero and N, we can find that: 

$ = fiN. (90) 

If the system has particles of various sorts numbered with i = 1, M and iVj is the number 
of particles of sort i, then in all the previous formulas the term fidN needs to be replaced 

M 

with ^2 fiidNi. We also have: 
i=i 

M 

$ = '%2ti i N i . (91) 

i=l 

In all of the above all the thermodynamic values, including the entropy, proved to be de- 
pendent on the values K[,...,K' k of the motion integrals K 1 ,...,K k . Now we're going to 
formulate the following 
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Equivalence principle. The entropy S(E, K[, K' k ) does not depend on K[, K' k (if 

all the other parameters of the system are constant). 

Here is how it can be substantiated. We assigned some specific value to the energy E, 
and for the sake of simplicity assume that n — 1. We normalize K' :— K[ in such a way as 
to ensure that K' changes from —1 to 1 and that at the entropy reaches its maximum. 

If S(K') depends on K' and our system is in thermodynamic equilibrium with its envi- 
ronment then it will be established such a value of K' that corresponds to maximum entropy, 
i.e. K' = 0. 

In a sufficiently small neighborhood of K' = 0, S(K') will look like this: 

S(K') = S -cK' 2 , (92) 

c > 0. Because our system is macroscopic and c must be comparable to So, we have to 
assume that c is very great. In this case the thermodynamic probability will look like this: 

W(K') = conste -0 *", (93) 

i.e. W(K') will have a pronounced peak at K' = 0, while all the other values of K' will 
be extremely improbable. This means that, in other words, that it doesn't really matter 
whether we're averaging using the micro-canonical distribution 

p(x) = c5(E - H{x)) (94) 

or using the generalized micro-canonical distribution: 

p(x) = c5(E - H(x))5(K' - K(x)). (95) 

Thus if S(E, K') depends on K' then we don't have to use the micro-canonical distribution 
with a fixed value of the integral K. 

Thus the most general case which we can come across is where the entropy S(H, K[, K' k ) 
depends only on the values of the integrals K 1 , K t , I — 1, k and does not at all depend 
on the values of the integrals i^z+i, K k . The integrals K-y, Ki are assigned values corre- 
sponding to maximum entropy and the mean for the generalized micro-canonical distribution, 
corresponding to the values of the integrals K x = K[,..Ki = K' t , K i+1 = K' l+1 ,...,K k = K' k 
will equal the meant for the distribution: 

k 

c5(H-E) J] 8(Ki-Kl). (96) 
i=i+i 
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The observed values of the integrals K i+1 , K k can be arbitrary. 

Now let us consider an example from the condensed matter physics where the general 
case above is realized. The case in hand is the superfluid helium. 

Let us consider a system of N particles enclosed in a certain macroscopic volume V and 
obeying the Bose statistics. The Hamiltonian of such system takes the form: 



N 2 N N 



^ = E|^ + EE $ (^-^i)' ( 97 ) 

i=l i=l j=i+l 

where — Xj\) is the potential energy of the i-th and j'-th particles. 
Let us assume that 

v(p) = J <5>(\ x \)e ipx d 3 x. (98) 
Then, in secondary quantization representation the Hamiltonian will take the form: 

T:=H-^N = J2(^-m b P + 

p 

E Hpi + P2 - p[ - p>(pi - p[HXA>M- (") 



Pl>P2,p'i,P2 



Here A(x) is a the function equal to 1 when x = and equal to zero in all the other cases. 
b^, b p are particle creation-annihilation operators in a state with the momentum of p. 

With low interaction between the particles almost all of them are in condensed state. And 
as b , 6q ~ y/~N, in commuting relation we can neglect 1 and consider bo, b$ commutation 
variables. Let us denote the number of particles in the condensate by A^. 

Keeping only the terms that are quadratic by the creation-annihilation operators of the 
super-condensate particles in the Hamiltonian, we derive: 

+^J2^pM b - P + b -P b p)- ( 10 °) 

P^O 

This Hamiltonian can be diagonalized by the canonical Bogoliubov transformation [13]: 

_ bp — Apbt p 

^ : " Trap 

b+ - Apb_ p 
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where 



= [N^W^Z (102) 
\l V m Am 

In new variables the Hamiltonian of the system will take the form: 

h = h + J2e p ^Z p - (103) 

If the condensate moves at non-zero rate u, then the Hamiltonian of the systems takes the 
form: 

H = H + J2( E p ~ u P)&Zp- ( 104 ) 

This Hamiltonian has the evident motion integrals — the occupation numbers. Or rather 
let us separate the set of all non-zero momentums {p} to disjoint subsets Si, containing Gi 
elements, corresponding to the close values of the momentums. We assume that Gi — > oo 
when V — > oo, but, together with that, ^ — > when V — > oo. Then the motion integrals 
will be the occupation numbers iVj of the cells Si. Let us assume rii — 7^. The entropy of 
the system will take form: 

S = ^Gi[(l + ni) ln(l + m) - n t In m}. (105) 

i 

This value reaches the maximum under the fixed energy and number of particles on Bose 
distribution. 

On the other side, except the occupation numbers, Bose gas has one more first integral 
- the condensate motion velocity. Let us analyze the dependence of the thermodynamic 
values of our system on the value of this integral. The state of the liquid helium is convenient 
to be described via the condensate wave function S [14]. It satisfies the following differential 
equation: 

i^(x, t) = -^-V 2 H(x, t) + U E(x, t)\E(x, t)\ 2 - cE(x, t), (106) 
where Uo, c > 0. 

If the system is kept at a constant temperature and its state is described by the number 
of Ai, A;, then in equilibrium the free energy F(T, X 1: A;) must reach the minimum. In 
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our case the condensate wave function Z(x,t) performs the role of these parameters. But 
the condensate wave function can vary with time and in meaning, Z(x, t), t) represent 

canonically conjugate variables. That is why the requirement of free energy minimality 
should be replaced by the requirement of stationary action: 

*2 

A[Z, S+] =J[iJ d'xE^E - F[Z, ~+]]dt. (107) 

As F is defined accurate to constant, we will find: 

F[Z,Z + ] = J d 3 x{-^Z+(x)V 2 Z(x) + ^\Z(x)\*-c\Z(x)\ 2 }. (108) 

To show that the thermodynamic characteristics do not depend on the velocity of the con- 
densate motion it is sufficient to show that the action is invariant under Galilean 
transformation: 

Z(x, t) ^ e i(2 ^ t+mra) S(x - vt, t), (109) 

v is relative motion velocity. But the invariance of the action A[Z, Z + ] in relation to the 
Galilean transformation is easy to find out by direct calculation. 

So, the action A[E, S + ] does not depend on the condensate flow velocity and this fact 
provides the stability of the condensate flow. 

We needed this example to illustrate the stability of the generalized microcanonical dis- 
tributions. Besides we obtained the answer to the question of why not all the first integrals 
residing in the system are observed. Only those integrals are observed that have the en- 
tropy unchanged while the values of these integrals are changing (at constant values of other 
integrals) . 

In the next section when we analyze the Ling cell, it will be essential that some of the 
first integrals of the Ling cell satisfy the equivalence principle. In addition to that, these first 
integrals have never been observed directly the way that the velocity of the Bose condensate 
through the capillary can be observed. These first integrals could only be observed indirectly 
by analyzing physiological phenomena testifying their presence. For this reason there would 
be a lot of doubts about or theory if no inorganic systems could be found for which there 
exist non-trivial first integrals satisfying the equivalence principle and which can therefore 
be observed. The example cited above allows us to dispense with such doubts. 
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Let us discuss the equivalence principle more detailed. For the sake of simplicity we will 
assume that the entropy of the system S(E, K') depends on the value of K' of one integral 
K, however, the reasoning given below is directly generalized for the case when the entropy 
depends on the values of the two and more first integrals in involution. 

According to the equivalence principle the entropy S(E, K') does not depend on K' . But 
it would be more precise to say that the entropy reaches its maximum on the certain interval 
of the values of the integral K' [a(X), /3(A)] and steeply decreases out of this interval. Here A 
- the parameter on which the Hamiltonian of the system depends. Here the energy of the 
system depending on the parameter A is chosen the way that at any value of A the entropies 
of the system were always the same. 

At the adiabatic infinitesimal change of the parameter A, A A' = A + SX the value 
of the integral K' is transformed by the formula K' h-» K' + §j<5A. Let us show that at this 
process the plateau [a(X ), f3(X )] transforms into the plateau [a(X'), /3(X')]. 

Our system is in equilibrium with the thermostat and let us assume that K' e e [a(Ao), P(Xq)] 
any equilibrium value of the integral K' (A = Ao). Let us show that when the system's energy 
is fixed at E and K 1 e [a(A ), /3(A )] the temperature of the system does not depend on K' . 

Let us assume that T — the temperature corresponding to the equilibrium value of K' e 
of the integral K, and equal to the temperature of the thermostat. Let us assume that 
K[ - - any value of the integral K from the interval [a(Xo), /3(Xq)]. Let us assume that 
T' — The temperature corresponding to it, and T' ^ T. Let us surround our system by 
an adiabatic membrane and change K' e h-> K[. When this happens neither the system 
energy + thermostat nor the system entropy + thermostat are not changing. Now let us 
make our system contact the thermostat again which will lead to the equalization of the 
temperatures of the system and the thermostat. At this stage we assume that system + 
thermostat is an isolated system, particularly, its full energy is constant. As at the process 
of temperature equalization a self-induced heat flow takes place, the entropy of the system 
and the thermostat increases. That means, at the value of the integral K K' = K' e for fixed 
full value of the common energy of the system and the thermostat, the entropy of the system 
and the thermostat is not maximal which contradicts the fact that K' e — the equilibrium 
value of the integral K. 

The same way we will demonstrate that at the fixed energy of the system the pressure 
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Pi, corresponding to the value of the integral K[ Is equal to the ambient pressure P. The 
pressure is defined by formula P = (§f )s- Let us assume, for definiteness, that P 1 > P. As 
it was demonstrated above the temperature corresponding to the value of the integral K[ is 
equal to the ambient temperature T. Let us make the system "expand" isothermally, doing 
work against the external forces. Meanwhile the system pressure Pi will be decreasing by 
virtue of the known thermodynamic inequation, fair in our case as well, that states that the 
module of the isothermal compressibility is positive. We will allow the system to expand 
quasistatically until its pressure is equal to the pressure of the thermostat. At that the 
system + thermostat will do a certain positive work, i.e. the full energy of the system 
+ thermostat will decrease by this value. But as the process of the system expansion is 
quasistatic and therefore reversible, the full entropy of the system + thermostat will not 
change. Therefore we can conclude that at the value of the integral K' e , and for fixed full 
entropy of the system and the thermostat, the energy minimum is not reached which means 
the value of the integral K' e is not equilibrium. 

So, we demonstrated that for fixed energy of the system E, The values of the pressure 
and temperature are constant on the whole interval [a(Ao), /?(Ao)] of the values of the integral 
K'. 

Let us now adiabatically change the A, A A' = A + 5X. At this change the entropy of 
the system will not change and the energy of the system will increase by the value PSX, i.e. 
the value that does not depend on K' . In other words, at the transformation K' i— > K' + ^-SX 
the points of the plateau will pass again into points of the plateau but corresponding to 
the changed parameter A. When making the inverse adiabatic change of the parameter A 
A' i— > A = X' — 5X, we see that the formula K' (->■ K' + ^-SX sets a one-to-one correspondence 
between the points of the plateau [a(A ), /3(A )] and [at(X'), f3(\')]. 

We remind that we proceed from the assumption that the system is in equilibrium 
with the environment but is not in equilibrium with itself, being in a stationary non- 
equilibrium state. Let us describe a possible way for the system to change from stationary 
non-equilibrium to equilibrium. Let us use the f)7Up for the distribution function. This ex- 
pression uses k independent first integrals in involution K\,...,Kk. However that does not 
mean that there no other integrals in the system except Ki, Kk, being in involution with 
them. For example in the formula for the ordinary microcanonical distribution only the 
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energy function takes part but, as it was discussed in the section 1, even the most realistic 
systems of the statistical mechanics are nonergodic at the thermodinamyc limit, i.e. there 
are other first integrals except the energy. We will call first integrals that are included in 
formula ( 1701 active. A possible way of changing to equilibrium is that the more and more 
first integrals from the list K x , Ky. cease being active. 

Let us assume that in this system the active integrals were K x , ...,Kk and the integral 
K k ceased being active. Let us use the symbol (SS) e to denote the change of the entropy at 
this process (at constant energy). Let us prove the following important proposition: 

Proposition. (6S)e does not change at quasistatic processes, that is to say 

(d(5S) E ) qs = 0. (110) 

Proof. Let us assume that the active integrals are K x ,...,Kk. Then the generalized 
microcanonical distribution corresponding to these integrals will take the form: 

k 

w k (E, K[,..., K' k , x) = cA(H(x) -E)H A(K t (x) - K[). (Ill) 

8=1 

Let us note that uik(E, K[, K' k , x) is in proportion to the characteristic function of a certain 
set. If the integral K k becomes inactive then the generalized microcanonical distribution oi 
Wk-i(E, K[, K' k _ 1: x) turns out to be composed of a certain number M of microcanonical 
distributions Wk(E, K[, K' k , x), the entropies of which are the same. It is evident that 
(5S) E = hiM. 

Let us consider the adiabatic change of the parameter A i— > X' = Xq + SX. We can assume 
that Vi = 1 , . . . , k 



(?r) , >=ao=o. (112) 

OX K[,...,K' k 

In this formula the low indices K[,...,K' k mean that the averaging corresponding to 
the vinculum is taken over the microcanonical distribution Wk{E, K[, K' k , x). This can be 
reached by replacing Wi = l,...,kKi h-> iiQ — (A— A )/i(i^i, K k , H), where fi{K x , K k , H) 
is the suitable functions of the first integrals and Hamiltonian. 

We will suppose that at this adiabatic process the point in the phase space that represents 
the state of our system, constantly moves in such a way that the values K x , ...,K k are the 
motion integrals. Meanwhile the microcanonical distribution Wk(E, K[, K' k , x) transforms 
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into the microcanonical distribution Wk(E', K[, K' k , x). But, because of the above, the 
value 8E := E' — E will not depend on the value K' k and will be the same if we calculat 
it for the distribution Wk-i(E, K[, K' k _ x , x). In the same adiabatic process the micro- 
canonical distribution Wk-i(E, K[, K k _ 1: x), corresponding to the value of the parameter 
A = Ao, will transform into the distribution w^i(E',K' 1 , ...,K' k _ 1 ,x), and by virtue of the 
fact that 5E can be calculated over the distribution Wk-i(E, K[, K' k _ x , x), the entropies 
corresponding to the distribution Wk-i(E, K[, K'^^x) with A = Ao and the distribution 
Wk-i(E', K[, ... : K' k _ 1: x) with A = A' will be the same. But at the adiabatic process (such 
as the one described in the beginning of the proof) the evolution of the distributions simply 
resolves into the canonical substitution of the arguments of these distributions. That is why 
with A = A', the distribution Wk-i(E', K[, K' k _^ x) will turn out to be composed from the 
same number of M generalized microcanonical distributions Wk(Ey, K[, ...,K' k ,x). That is 
to say M = const. That is to say ((SS) E ) qs = InM = const, which was to be proved. 

In case if there are too many of first integrals and if the process of decrease of the number 
of active integrals can be considered continuous, it is reasonable to introduce the continuous 
parameter s, characterizing the number of active integrals the way that s G [0, 1], the value 
s = corresponds to the case when no integrals are active and the decrease of the number 
of active integrals corresponds to the decrease of the parameter s. Therefore we have: 

(113) 

The formula ^ < is correct when it is assumed that the system energy and volume are 
kept constant. Let us assume, however, that the system temperature and volume are kept 
constant. Then: 

dS _ dS dS dE 

- { i; )Ey ~ { Ts )Ty ~ { dE )vA i; )v ' T ~ 

dS 1 ,dE 

1 AST-El 
~ T ( Is hv - (U4) 



That is to say: 

(^)t,v<0. (115) 

The same way let us assume that the system temperature and pressure are kept constant. 
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Then: 



dF. .dF. dF. ,dV 




dF, „,dV 



(116) 



So, in case when the temperature and pressure are kept constant: 




(117) 



The statement, that the number of independent first integrals in involution for the system 
must be very large can be established, for example, by the following way. Let us decom- 
pose the system & into disjoint union of enough large number N of small but macroscopic 
subsystems ©j, % = 1, ...,N. Let H iy i = 1, ...,N be a Hamiltonian of subsystems (3j and 
Hij, i,j = 1,...,N, i ^ j be the total interaction Hamiltonian of subsystems &i and &j. 

N ' N 

Then Hi is proportional to the volume of system & and ^2 Hij * s proportional to 

i=l l<i<j<N 

the total area of boundaries between different subsystems ©j. But each subsystem &i is 
macroscopic. So we can neglect by the interaction Hamiltonian of different subsystems in 
total hamiltonian H of system &. Therefore we can write that 



But each subsystem &i is macroscopic. So, according to the nonergodic theorem for each 
subsystem (3, there exists at least one its first integral Ki. It is easy to see that for different 
i = 1,...,N integrals Ki are independent and in involution. It follows from fl 11 8 [) that 
integrals Ki are the first integrals of whole system & at the same time. 

5 E.S. Bauer's stable non-equilibrium principle and 
the resting state of G. Ling. 

For the sake of further analysis we will consider the contribution that E.S. Bauer and G. 
Ling made to the issue at hand. According to Bauer [4] the substance of the living cell can 
be in one of two states: the stable-non-equilibrium (resting) and thermodynamic equilibrium 



N 




(118) 



i=l 
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(active). The work The work produced by living system a living system is done when the 
system passes from stable non-equilibrium to equilibrium. Bauer used this principle (the 
two basic states and the change between them) as the basis for his "theoretical biology" and 
deductively derived the basic properties of such biological processes as assimilation, growth, 
excitability, adaptability and reproduction. That is what Bauer himself says [4, p. 143] 
about this principle: 

"Nonliving systems are ever in equilibrium and due to their free energy they do work 
against equilibrium required by the laws of physics and chemistry in the existing ambient con- 
ditions." Further Bauer expands on it: "We will call this principle as stable none-equilibrium 
principle of the living systems. This name clearly expresses the meaning of the principle and 
the unique feature of living systems from the point of view of thermodynamics. The same 
way as a system in stable equilibrium when disturbed returns back to it, a living system 
in stable non-equilibrium also keeps coming back to it. Our principle also characterises the 
unique feature of living systems as we do not know any nonliving systems that would be in 
stable non-equilibrium." 

The ground rule of Ling's Association-Induction Hypothesis [1,2,3] is that the cell as a 
quasi-solid body is a system with a non-maximal entropy in a resting state. We proceed from 
the approach based on the fact that the resting state of the living cell, according to both 
Ling and Bauer, is a stationary non-equilibrium state, the possibility of existence of which 
was proven by one of us for a large class of realistic systems in statistical mechanics [7]. Here 
we simply apply it to the living cell. The appropriateness of the proposed approach to the 
living will be verified by comparing the results of the theoretical analysis and the available 
experimental data. 

Based on the assumption that the state of the living cell is stationary non-equilibrium 
we will give a thermodynamic description of the following phenomena: 

1) When the cell is excited and when it dies heat is generated. 

2) When the cell is excited and when it dies the cells' size changes (mostly decreases). 

3) When the cell is excited and when it dies , the cell's key protein molecules fold. 

4) When the cell is excited and when it dies, efflux of potassium ions from the cell takes 
place. 

Let us begin with the explanation of the first phenomenon. Let us assume that the 
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substance inside the cell is in stationary non-equilibrium corresponding to a number of active 
first integrals in involution and the process of excitation and death manifests through some 
of these first integrals becoming inactive. Let's assume that just like in the section 3 that 
the number of the active first integrals is characterised by a continuous s G [0,1], Where 
s = corresponds to the case, when neither of first integrals are active and the value of 
s = 1 to the case when the maximum number of the first integrals are active. The question 
why there must be so many nontrivial integrals that their number can be described by a 
continuous parameter, can be answered by the fact that the muscle fibres can be contracted 
to different degrees and is able to algebraically summate the nerve impulses that enter it. 
It will be clear when we consider the changes in the size of the Ling cell when it's activated 
and when it dies. We can choose s in such a way that s is proportionate to the number of 
active first integrals. The system's volume is supposed to be constant. 
Let us make an infinitesimal change in the parameter s: 

s ' = s-6s, (119) 

where 5 > 0, 5 s is are infinitely small. 

It has been established that the knowledge of k of the first integrals in involution allows 
decreasing the number of degrees of freedom by k units (system deflation according to Whit- 
taker). For each s we have const s of active first integrals and the entropy that corresponds 
to the parameter s is the entropy of the reduced system, corresponding to these constant 
s first integrals. If s s i— > s' = s — 5s, then the number of degrees of freedom of reduced 
system increases by consttfs units, and we may consider that these new 'turned-on' degrees of 
freedom corresponding to the new system 5&, that is in thermodynamic equilibrium with the 
initial system &. We will change the parameter s at a fixed temperature T. If s (->■ s' — s — 5s 
the system's energy changes: 

E' = E+(5E) T , (120) 

and entropy: 

S ^ S' = S +{5S) T . (121) 

We can equate (5E) T to the energy of the system 5&, and (5S)t to the entropy of the 
same system. As the system 5& is in thermodynamic equilibrium with the system &, it is 
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essential that the temperature of the system S& should also be equal to T, in other words, 
it is necessary that: 

1_ _ d(SS) T 
f ~ d(8E) T 

But this equality can also be derived directly. We have: 
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dS _ 1 

Ie~ f'' 

dS + d(5S) T 1 



(123) 



dE+d{8E) T T' 

From these two equalities we derive (11221) by using the elementary transformation. 
But in order for 5& to be in equilibrium with & the equation (11221 is not enough. It is 
necessary that 

gjgk < (124) 

This inequality results from the requirement for maximum entropy in the system (3 + 5& 
under condition that the energy remains constant, if we take into account the fact that the 
number of degrees of freedom of the system 5& is far less than the number of degrees of 
freedom of the system &. The detailed discussion of this inequality is found in Appendix 1. 
This inequality leads to 

d i<0. 



d(8E)x T 

' ,IT <0, (125) 



T 2 d(SE) T 
or 

:= d -^k > o. (126) 

We will make the following assumption about the value of (5E)t- There is a temperature To, 
that makes the temperature T > TO incompatible with the existence of life in the cell. This 
means that if the temperature T — > Tq the interval A from the formula f!70|) becomes more 
and more comparable to the range spaces of the integrals Ki, ...K^ values, and the generalized 
micro-canonical distribution degenerates into a usual micro-canonical distribution. Then, if 
T>T 

(5S) E = const, (127) 
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) i.e. if T > T the dependence of T on E for & will be the same as the one for & + 5&, 
i.e. (5E) T = 0. This boundary condition is discussed in more detail in the Appendix 2. We 
have: 

To 



(SE)t = («£)* - / ^<£T< 



T 

T 



= (5E) T() - J 5c v (T')dT' = 



T 
T 



~y Sc v (T')dT' < 0. 



(128) 



T 



I.e. (5E) T < 0, and that's what we set out to prove initially. 

Let us proceed now to the explanation of why the Ling cell changes size when it is excited 
and when it dies. Let us examine a resting cell in equilibrium with the environment at a 
constant temperature T and pressure P, and show that when the cell is excited and when it 
dies its size decreases. So, we assume that as above s h-> s' = s — 5s, 5s > 0, and that (SS)e 
means the change in the entropy of the cell if its energy is constant. 

In section 3 we demonstrated that (5S) e = const in an adiabatic process and that is why 
(5S) e is a function of the entropy only, that is to say: 

(5S) E = 5f(S), (129) 

for some function 5f of one variable. Let us show now that 5f is a decreasing function. We 
have: 

(5S)e = (5S) t -(5E) t —, 
(5S) e =(5S) t -^(5E) t . (130) 

Because of the thermodynamic inequality cy > ((cy -the thermal capacity of the cell 
at constant volume) the energy, and as a consequence, the entropy are strictly increasing 
functions of the temperature. That is why it is enough to show that -^p((5S) E ) < 0. But 

^((SS) E ) = ^((SS) T ) - ~{{SE)t) + ^(SE)t = 

= ^(SE)t- (131) 
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And by virtue of this, as it was already proved, (5E) T < we find: J^((5S)e) < 0, i.e. Sf 
is a decreasing function of its argument. 

Now let us define the minor additions that appear in the thermodynamic potentials under 
s (->■ s' — s — Ss. We have: 

(SS) E = (SS) T -(SE) T —, 
(SS) E =(SS) T -^(SE) T . (132) 

I.e. under s h-> s' = s — Ss, 



E(V, S) ^ E(V, S) - T8f(S). (133) 

According to the minor additions theorem for free energy F and the thermodynamic potential 
$ we have: 

F(V,T)^F(V,T)-T5f(S), 

$(P,T) ^$(P,T) -T5f(S). (134) 



As 



v =( — Qp — )t, (135) 



Then change in the cell size with s \— > s' = s — 5s takes the form: 



r)V r)^ 

5V = -T(5f)>(S)(— ) T (— ) T . (136) 



As (5/)' < by virtue of what was proven above and (§^)t < (thermodynamic inequality, 
which is always true), the sign 5V is opposite to the sign (§^f)t- 

Let us set the sign (§^f)t- But because of a known thermodynamic relation 

,dS. ,dV. f dP, , i0 _ 
(^)t = -igfM W )T- (137) 

That is why the statement that (f^)r > is therefore equivalent to the statement that 
the cell expands when heated (at constant pressure). This behavior is demonstrated by most 
solid bodies when heated. If we accept that this is true (the cell expands when heated at 
constant pressure) then we find 5V < 0, which was to be proved. We thus conclude that the 
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cell's size decrease. We consider our description to be a general thermodynamic prerequisite 
that explains the phenomenon of the Ling cell contraction, a prototype of muscle contraction. 

Before we consider the aggregation of key proteins in the Ling's model, let us go back to 
the question of the transformation of the internal energy of the cell when it dies. We have 
demonstrated above that death occurs then 

{8E) Ty < 0. (138) 

If we use the symbol 5Q to denote the amount of heat that entered the cell, we will find: 

dS 

5Q = (5E) Ty + T(—) T 5V. (139) 
But we already know that SV and have opposite signs, i.e. 

5Q < 0. (140) 

In other words, the death and excitation of the cell are exothermic reactions. 

Let us consider now the thermodynamic prerequisites of the aggregation of proteins that 
are present in a resting Ling's cell in an expanded state. 

Our analysis draws on the fact that the energy of an unfolded protein molecule is greater 
than the energy of a folded protein molecule and is based on the protein model in which the 
configuration of the protein molecule can be characterized by a real parameter x e [0,1], 
describing the degree to which the molecule has unfolded. We assume that the energy of 
the protein molecule is a function only of x and it is an increasing function. We assume also 
that x = corresponds to a completely folded molecule whose energy is at a minimum; and 
x — 1 corresponds to a completely unfolded molecule whose energy is at its maximum. 

Let us use the following supporting procedure. We will identify a single protein molecule 
in the cell's protoplasm and call it X, and we'll ULS6 X clS du parameter characterizing the 
degree to which our single molecule is unfolded. If we call the system's temperature T, then 
all the thermodynamic equations derived in the previous section can be applied to X, more 
particularly we have the equation: 

dE = TdS - hdx, (141) 

for a certain function h(E, S) of the energy and entropy X. The value h(E, S) is a generalised 
force that is thermodynamically conjugated to x. When the external forces that affect our 

43 



protein molecule, are absent, h(E, S) = 0. The relation (I14ip is equivalent to the equation 

dE = TdS - PdV, (142) 

And that is why we can more or less directly apply to case the conclusions for the cell 
interacting with the environment, changing x ^ V, P ^ h. So, Let's assume again that 
s i— > s' = s — 5s, then 



(5x) T , h = -T(5f)'(S)(—) T , x . (143) 

We see that the sign of 5x is the same as the sign of (§f )t- by virtue of the known 
thermodynamic relation: 

iff* --<!)». ("*) 

So, the sign (|f )t = — (ff)ft wm on ^y be negative if the protein molecules unfold and the 
cell temperature increases. So we have proved that the protein molecules fold when the 
protoplasm is activated and when the cell dies only if they unfold when the cell is heated. 
It is clear that this statement is true only for the protein model used here. 

But the property of the protein molecule to unfold when heated appears to be physically 
evident. This assumption seems to be reasonable as the protein molecules are supposed to 
change to states with more energy when heated, and the energy is an increasing function of 
the parameter x. That is why the protein molecule folds under s s h- > s' = s — 5s, which was 
to be proved. 

The question of how we apply this general method we used for the analysis of more 
realistic protein models needs further research. 

Let us consider now the question of why the Ling cell's nonergodicity includes ion flows 
between the cell and the environment when the physiological state of the Long cell is changed. 
When analyzing this phenomenon we will restrict ourselves to potassium ions that play an 
important role in the physiological processes of the cell. We will try to answer the question of 
why this positive ions are released by excited or dead cells. The answer to this question will 
explain another phenomenon well-known from the classical cytology: why the damaged area 
in a cell always has the negative charge to the adjoining intact areas that are not affected 
by excitation or injury. 



44 



It is clear that our approach must be improved in order for us to be able to consider a 
broader spectrum of ion characteristics, such as: ion chemical nature, their physical proper- 
ties, their concentration in solution, etc. 

Let us assume that N means the number of potassium ions in the cell <£. We can 
regard iV as a parameter, and it will be now clear why. Our analysis is carried out in the 
thermodynamic limit and from this point of view the cell can be considered very big. Let 
us attach a very small vessel (3, to the cell which will be separated from it a semipermeable 
membrane that is permeable only for potassium ions. We will assume that in & we have 
only potassium ions and that & is so small that the entropy of the ions inside it is equal 
to zero. As a model for & we can take a very narrow potential well, the quantum energy 
levels of which are not degenerate and the distance between adjoining energy levels is much 
greater than the temperature T. Then potassium ions in & will be mostly in a basic energy 
state and the thermodynamic probability W — 1, and S e = \nW = 0. Let us also introduce 
an external electric field the potential of which ip everywhere except & is equal to zero, 
and inside & is constant and equal to ipo. Adiabatically changing tpo, we can pump the 
potassium ions from & to the cell and back. Meanwhile the entropy of the whole system 
will be constant and equal to the entropy of the cell S. By virtue of the fact that <p is 
a parameter in the Hamiltonian of the whole system, ((SS)E) qs — const during adiabatic 
change ipo. Here E means the energy of the entire system. Let us assume that E<r means 
the energy of the cell. All the motion integrals that are related to the entire system & + € 
are equally related to the cell <£ because the potassium ions in & are immovable, and their 
coordinates and momentums can be removed from the list of the dynamic variables. But 

(SS)e = (SS) E£ , N + ^-8E € + —8N. (145) 

But 8 En = —fo8N, as the full energy of the cell and of the vessel is constant. Considering 
that in equilibrium: 

-3E^ m+ m SN=0 - (146 > 

we find: 



(8S) E = (8S)e € ,n- (147) 
That is why during an adiabatic change in the parameter <p or, of the number of particles 
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in the cell N 



((SS) Ec ,N)qs = const, (148) 

Which means that N can be considered as a parameter. 

Let us proceed now to the question of the amount of potassium ions released by the cell 
during its activation and death. We assume that the cell is immersed in a potassium ion 
solution and the chemical potential of the potassium ions is constant and equal to fi. The 
following thermodynamic relation will hold true: 

dE = TdS + fidN (149) 

and all its corollaries. This equation similar to the thermodynamic equation 

dE = TdS - PdV, (150) 

if we equate N ■<->■ V, /i ■<-)■ —P. Therefore, to prove that the potassium ion yield from the 
cell takes place when s H- s' = s — 5s all we have to show is that 

(m )t>0 - (151) 

To find (§^)t we will note that the concentration of potassium ions in the cell is low and 
that is why we can use the strong electrolyte theory (see [13]). This theory suggests that 
in a low concentration of potassium ions the cell's entropy (at fixed temperature) S = 
S + S K+ + Nip(P,T), where S'o - the cell entropy at zero ion concentration, S K+ - the 
entropy of only the potassium ions, calculated for a perfect gas, and <p(P, T) is a function 
of the cell temperature and pressure. That is why for the entropy of the potassium ions we 
can use the standard formula (see [13]): 

S K+ = N\n— -Nf'(T), (152) 

where 

/(T) = -Tln((^) 3 / 2 ). (153) 

Here m is the mass of a potassium ion. It is clear that if the potassium ion solution in the 
cell is thin enough (ffjOr > 0. 

46 



So, 5N < 0, which was to be proved. 

For simplicity of the analysis we assumed above that there are only potassium ions in the 
cell. But in reality there are different kinds of ions and all the derived formulas should be 
true for each kind of ions. Particularly, during the activation and death of the cell all kinds 
of ions should come out of the cell if their concentration is low enough. But it is known 
that some ions, for example sodium ions, are on the contrary, absorbed by the cell during 
activation or death (according to Ling's model). . However, the cell becomes negatively 
charged during activation, i.e. our explanation leads to the correct answer "on average". 
This gives us hope that our explanation is genuinely true and can be improved in a proper 
way. 

To conclude of this section let us note that we considered the resting state of the cell to 
be stationary non-equilibrium and the state corresponding to the excitation and death to be 
equilibrium. We considered these states as given and never asked the question why there are 
transitions between these two states. One of the possible hypotheses is that such transitions 
appear due to changes in the conditions on the cell's boundaries. That is why in section 3 
we discussed in detail the role that the conditions on the system's boundaries play in the 
establishment of thermodynamic equilibrium inside the cell. The question of how true this 
hypothesis is will be the subject of further research. 

6 Conclusion. 

We used the results of [7] to analyse the Ling cell in this paper, assuming that this kind of 
cell is a nonergodic system. 

Using the property of nonergodicity as a starting point for our research, we first discussed 
the nonergodicity of statistical mechanics systems, we then developed the thermodynamics of 
stationary non-equilibrium states assuming the existence of several nontrivial first integrals 
that commutate between them. These thermodynamics are applicable to any system in 
statistical mechanics, including the Ling cell. 

The analysis of the Ling cell properly demonstrated that using the approach based on our 
thermodynamics of stationary non-equilibrium we managed to derive such realistic properties 
of the excited Ling's cell as heat generation, decrease of the cell size, the folding of the 
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proteins that are unfolded when the cell is resting and the releasing of potassium ions. The 
applicability of the proposed approach to the living cell model, such as the Ling cell, testifies 
to the applicability and adequacy of our analysis for studying the living state of matter. 

We are very grateful to P. Agutter, A.V. Koshelkin, Yu. E. Lozovik, A.V. Zayakin and 
E.N. Telshevskiy for valuable critical comments on this article and very useful discussions. 



7 Appendix 1. Discussion of the inequality (11241) . 



In this section we will try to justify the inequality ( I124p . 

Let us see the Hamiltonian dynamical system with n + k degrees of freedom (k <^ n) 
and k independent first integrals in involution Ki,...,K k . Locally we can always choose 
the canonical coordinates (pi, qi),...,(p n+ k, q n +k) in such a way that Ki, ...,K n will only be 
the functions of (p n +i, Qn+i),---,(Pn+k, Qn+k)- First, let us assume that the coordinates that 
measure up this property can be chosen globally. So, the phase space M of our system 
represents as a direct product of the phase spaces M = Mi x M 2 , where Mi corresponds to 
the coordinates (pi, qi),...,(p n , q n ), and M 2 — to the coordinates (p n+ i,q n+ i),...,(p n+k ,q n+k ). 
Let us denote through dT the element of the phase volume on M, and through dT 1 and 
dT 2 the elements of the phase volumes on Mi and M 2 , correspondingly. Let us assume that 
H(x, y) — the Hamiltonian of our system, x G Mi, y e M 2 . 

Let us assume that S(E) is the entropy calculated for the microcanonical distribution and 
Si(E, K[, K' k ) is the entropy calculated for the generalized microcanonical distribution 
corresponding to the K[, K' k of the first integrals Ki, K k . Let us recall that 

Si(E, K[, K' k ) = In I dT5{H - E) JJ 5{K, - K[). (154) 

J i=l 

We assume that Si(E, K[, K' n ) measures up the equivalence principle (see section 4). 

Locally on M 2 we can build the function tpi, <pk> an d do it in such a way that the set of 
functions (ifi, (pi),...,(K k , ip k ) will represent the set of canonical coordinates, i.e. the Poisson 
bracket in these coordinates will have the standard form. In other words, we can build a 
partition M 2 in the area Vi, i = 1,2, ... with a piecewise-smooth boundary and in each of 
them we can choose the functions <p\, <p\, that satisfy the requirement mentioned above. 
Let us assume that £ — a certain joint surface of the level of the integrals Ki, K k , and 
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x G E. If (without losing generality) we assume that E is connected, then the Hamiltonian 
phase flows generated by Ki, ...,Kk transitively acts on E. And as these phase flows keep 
the phase volume on M and the Hamiltonian, then 

S 1 (E,K[,...,K' k ) = \n J dT l x S(H(x,y)-E) + c, (155) 

where x G Mi, and c is the constant equal to: 

c = ln E J M-M- (ise) 

The constant c does not depend on E and can be discarded. 

Let us now turn to the calculation of (5S)t- If the system M = Mi x M 2 is described 
by the Gibbs distribution then the probability to find the system M at the point (x, y), x G 
Mi, y G M 2 takes the form: 

1 H(x,y) 

where 

Z (T):=/ W ^. (158) 
The probability to find the system M 2 at the point y is given then by the expression: 

w 2 (y)= I dTlw 12 (x,y). (159) 



Let us introduce the conditional probability Wi\2{x\y) to find a system Mi at the point 
x at condition that the system M 2 is at the point y by the formula: 

wi 2 (x, y) = wi\ 2 (x\y)w 2 (y) . (160) 

Von Neumann entropy of the system M takes the form: 

S{E{T)) = - J dT l x dT 2 y w 12 {x,y)\nwi 2 {x,y). (161) 

Using (11601) S(E(T)) we can transform to the form: 

S(E(T)) = (Si(E(T),y)) M2 + S 2 (T), (162) 

where 

Si(E(T),y) = - J dT l x w A2 (x\y)hiwi\ 2 (x\y), (163) 
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(164) 



Note that only averaged values can be measured at experiment. We find: 



(5S) T = S 2 (T). 



(165) 



At that we meant the formula ( 1155j) and the fact that within the thermodynamic limit the 
descriptions by the microcanonical and canonical assemblies coincide. At the same time we 
also consider, of course, that the constant c from the formula HI 56[) does not depend on 
K[,...,K' k . When we will proceed to the discussion of the general case that will not require 
the presentability of M in the form of the direct product Mi and M 2 we will see that this 
assumption is justified. 

Let us try to describe the dynamics of the system M 2 , using the fact that the number of 
degrees of freedom M 2 n is much greater that the number of degrees of freedom M\ k. The 
Hamiltonian canonical equations for M 2 take the form: 



As the number of degrees of freedom M 2 is much less that the number of degrees of freedom 
Mi we can average the right parts of the last equation by Wi\2(x,y). As a result, skipping 
some quite trivial calculations we will find: 



Pi = 



Qi = 



dH(p 1 ,q 1 , ...,p n+k ,q n+k ) 
dqi 

dH(p 1 ,q 1 , ...,p n +k,Qn+k) 
dpi 

% — n + 1, n + k. 



(166) 



Pi = 



Qi = 



9F(p n+ i, q n+ i, ...,p n+ k, q n +k 

\T) 

dqi 

9F(p n+ i, q n +l, ■■■,Pn+k, 1n+k\T) 

dpi 

% = n + 1, n + k, 



(167) 



where 



\T) 



TlnZ 1 (y\T) 




V — yPn+li 9n+l) ■■■■.Pn+ki Qn+k)- 



(168) 
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So, this way the dynamics of the system M 2 will be described by the canonical Hamilto- 
nian equations with the Hamiltonian F(y\T). 

Let us note that this way the defined Hamiltonian is not defined uniquely but accurate to 
the temperature arbitrary function f(T). Let us define the new Hamiltonian of the system 
M 2 

H 2 (y\T) = F(y\T) + f(T), (169) 
in such a way that makes true the following: 

(6E) T = (H 2 (y\T)) M2 . (170) 

Let us show that in this case 

(-^H 2 (y\T)} M2 = 0. (171) 



For this we will proceed from the identical equation: 

(172) 



But 



where 



(SS) T = S 2 (E(T)) = - j dT 2 y w 2 (y)\nw 2 (y). 

dT 2 y e-^. (174) 



We have: 



P^EiT)) = _(A(^I))A(A^ )>Mi , (175 ) 
where we assumed by definition: 

&g(y) ■= g(y) - (f(y)W (176) 

Making the similar calculations, we find: 

^(H2(y\T)) = (±H 2 (y\T)) ~ (A(H 2 (y\T))A(^^^)) M2 . (177) 
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Comparing (j!75p and (11771) we will find: 

(-^H 2 (y\T)) M2 = 0, (178) 

which was to be proved. 

Let us describe now the effective Hamiltonian of the whole system Mi x M 2 . The dynamics 
of the system M\ can be described by the Hamiltonian H(x, y) where y G M 2 and depends 
on time, and x G M\. But because of Ki, ...,K^ are the motion integrals, the value H(x,y) 
at fixed x is the same for all y corresponding to the same value of the motion integrals. That 
is why the system Mi is described by the Hamiltonian H(x,y), where y is constant with 
time and is chosen the same way as in the formula (j!55p . 

Further the effective Hamiltonian of the system M 2 H 2 (y\T) depends on the temperature 
T of the system Mi, but temperature is an intensive parameter and instead of T we can 
use the energy of the system Mi, falling within one degree of freedom. Let us assume by 
definition: 

H' 2 (y\^) = H 2 (y\T), (179) 
n 

where E\ — the energy of the system Mi. 

As an effective Hamiltonian of the system M = Mi x M 2 the following expression will 
be used: 

H eff (y, z) = H(y, x) + H' 2 (z\^^-). (180) 

Let us recall that according to this formula x is chosen the same way as in (I155p . 

Let us show that within the limit n — > oo the Hamiltonian equation for H e ff will coincide 
with the Hamiltonian equations for the initial Hamiltonian if the point x is on the same joint 
surface of the level of the integrals K 1: K k , that z. 

Let us assume that (y, z) G M = Mi x M 2 . Let us analyze the derivative with time of 
the point z by virtue of the Hamiltonian equations for the Hamiltonian H e ff. This will be 
exactly the derivative by virtue of the Hamiltonian motion equations corresponding to the 
Hamiltonian H 2 (z\T) and the fact that this derivative coincide with the derivative with time 
by virtue of the Hamiltonian motion equations built by H(y, z) was shown by us above. Let 
us assume now that R(y) is the dynamic variable on Mi. Let us calculate its derivative with 
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time by virtue of H e ff. We have: 

R(y) = (H(y,x),R(y))(l+-H';(z,^^-)). (181) 

n n 

Here the symbol H' 2 \s, z) denotes the derivative H' 2 (s, z) by the first argument. Within the 
limit n — > oo the right part of (j!8ip evidently transforms into (H (y , x) , R(y)) , i.e. within 
this limit 

R(y) = (H(y,x),R(y)), (182) 

which was to be proved. 

Let us note as well that H e ff depends on the choice of the point x G M 2 . However, 
because of the equivalence principle, the entropic properties of the system M = Mi x M 2 
will not depend on this choice and these last ones are the only ones that are important for 
us. That is why from the point of view of the calculations of the entropic properties the 
Hamiltonian H e ff(y, z) is as good as the initial H(y, z) and we will further use only the first 
one. 

Let us note finally that the effective Hamiltonian with the form of (11801) is defined accurate 
to the arbitrary function / ( H ^ y,x ^ ) . I.e. the effective Hamiltonians corresponding to the 
different choices of such normalization induce the same motion equations. That is why all 
the thermodynamic characteristics calculated by H e ff corresponding different choices of the 
normalization must coincide. And this is really so. Actually at the old temperature the 
adding of f( H ^' x ^ ) to H e jj the arbitrary function can be considered by the corresponding 
change Hx(y, x). But at the same time with n — > oo the average value of the energy E\ of 
the system M\ will change by the value ~ 1, and that is why the entropy of the system M 2 
will change by the value ~ i i.e. it will not change at all within the thermodynamic limit. 

Let us use S(T) to denote the entropy of the system M, as the function of the temperature 
T, calculated using Gibbs distribution for the Hamiltonian H e ff. Let us use E\ to denote 
the average Hi(y,x), and E 2 to denote the average H 2 (y\T) by the measures on Mi and M 2 
respectively, induced by the same Gibbs distribution. Let us assume also E = E\ + E 2 . The 
direct calculation shows that ^| = ^. From the other point of view we normalized H 2 {y\T) 
in such a way that dS2 ^J^ 2 ^ = ^ would be true. From these two equalities we obtain 
^dE^ = T- I- e - m ^ ne thermodynamic equilibrium state the temperatures of the systems 



53 



Mi and M 2 coincide, and this equality is true precisely and not only in the thermodynamic 
limit k = const, n — > oo. 

If k is very large (being at the same time much less than n), then according to the main 
principles of the statistical mechanics, instead of the description of M 2 by means of Gibbs 
canonical distribution we can use the description by means of the microcanonical distribution 
and the function of the distribution for the whole system M, at which the energies of the 
systems M x and M 2 will be equal to E 1 and E 2 respectively, will take the form: 



To obtain the inequality that we need let us make the following conceptual experiment. 
Let us make the system M 2 contact the thermostat T, which is at the temperature T. As 
for the system M 1 , we will take a very large system M 3 that will be connected to the system 
Mi via a heat pump N, working by Carnot cycle. We will consider the working body of this 
pump to be so small that its thermal capacity can be neglected. Let us also assume that 
during the thermal contact with Mi or M 3 the temperature of the working body of the pump 
is equal to the temperature of Mi or M 3 (depending on which one has thermal contact with 
the working body of the pump), i.e. the heat pumping is done without the entropy increase. 
The working body of the pump is connected to a shaft by something like a connecting rod 
gear, this shaft can rotate and on the axis of which there is a spring P, possessing the energy 



Let us note that because of our assumptions about the effective Hamiltonian H e ff(y,z) 
the dynamics given by it on M 1 are exactly the dynamics given by H(y,x) and that is why 
the system Mi can be considered closed. Or more precisely as the system Mi is related 
to M 3 by means of the heat pump N the system that consists from Mi, M 3 , N, P can be 
considered closed. We will now assume that the elastic characteristics of the spring P are 
picked the way that the work that must be done to the working body N is exactly equal 
to the change of the spring energy. So this way we have a whole continuum of equilibrium 



5{H{y,x)-Ei)8{H' 2 {z\^)-E 2 ), 



(183) 



y G Mi, z G M 2 . The entropies of the systems M x and M 2 will be equal to 




(184) 



Ep. 
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states of the system consisting of M\, M 3 , N, P, which is parametrized by, for example, the 
angle of rotation of the pump shaft and at Si + 63 = const, E 1 + E 3 + Ep = const. 

When our whole system obtains the equilibrium then the temperature of the thermostat is 
equal to the temperature of Mi and the temperature of M3. From the condition of maximal 
entropy of the whole system that consists of Mi, M 2 , M 3 , T, N, P at fixed full energy 
considering the equalities: Si + S3 = const, E x + E 3 + E P = const we will find that in the 
equilibrium the following value 

S 2 {E 2 , Ei) + S T (E T ) (185) 

should reach the maximum at the additional condition 

E 2 + E T = const. (186) 

Here Et — the energy of the thermostat and St(Et) — the entropy of the thermostat. 

Instead of the dependence S 2 (E 2 , Ei) on Ei it appears more convenient to consider the 
dependence S 2 from the system temperature Mi X(Ei). Let us introduce a new function 

S'(E,\(Ei)):=S 2 (E,Ei). (187) 

The condition of the maximality of the entropy of the system M 2 and the thermostat 
lead to the equalities: 

dS'(E, A) 



dX 



\\=T, E=E(T) — 0, 



dS'(E, A) _ 1 

7^ \\=T,E=E(T) - (188J 

But these equalities can be obtained directly. The second equality is obtain from the ordinary 
formulas of the statistical mechanics if we assume that the system M 2 is described by the 
Gibbs canonical distribution with Hamiltonian H 2 (y, A). Let us set up the first equality. Let 
us adiabatically change the parameter A, A H- T + 5X. At such change of the parameter A 
the entropy of the system M 2 will not change. But by virtue of (11711) the energy M 2 will not 
change. I.e. the derivative of the entropy of the system M 2 by A at fixed energy is equal to 
zero, which is exactly the equality that we needed. 

But for the maximality of the entropy of the system M 2 and the thermostat it is also 
necessary that the matrix of the second derivatives S'(E, A) at the point E = E(T), X = T 
would be defined negative. 
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To prove now the inequality (11241) . we should prove that 

d 2 



dE 2 



S'{E,T{E)) < 0. (189) 



But, 



dE 2 v ' v " dE 2 dE OEdT 

And the last inequality is true by virtue of the negative determination of the matrix of the 
second derivatives S'(E, A) at the point E = E(T), X = T. Which was to be proved. 

In a general case when M is impossible to be conceived as a direct product M = Mi x M 2 
it can be shown that for some covering (M, H) of our Hamiltonian system (M, H) such 
factorization is possible. Deriving from the very beginning all the thermodynamic relations 
for (M, H), as it was done in section 3, we come again to the inequality (11241) . 

Or more exactly we want to say the following. In all our thermodynamic analysis the 
Hamiltonian H and the integrals Ki, ...,Kk depended on a certain parameter A ( volume), 
H(X), Ki(X), K k (X). Further the dynamic variable X(X), depending on the parameter A 
at A = will be denoted simply as X. 

We can build a covering Hamiltonian system (M, if) of the system (M,H), so that M 
would be conceived as a direct product M = M\ x M 2 and so that will measure up one more 
additional condition. Let us assume that n is a canonical projection of M on M and let us 
assume that H(X), Ki(X),...,K^(X) are the lift H(X), Ki(X),...,Kk(X) on M. The additional 
condition mentioned above is that the canonical coordinates on M2 are K\, K k and the 
conjugated variables (fx, (fk, and Wi = 1, k tpi runs all the real axis. 

We will now derive all the thermodynamic relations for M, H. But here we get the 
difficulty related to the fact that M2 is not compact. There is how we propose to overcome 
this difficulty. We propose to deal with unnormalized distributions w for which the following 
is true J w(x)dT x = 00, where dT x is the element of the phase volume on M. 

M 

Let us assume that Rl = {x e M\Wi = 1, k \<fii(x)\ < L}. To calculate the average of 
the dynamic variable / on M by w we propose to use the following formula: 

/ f(x)w(x)dt x 

(/> = lim ^— = . (191) 

L ^°° J w(x)dT x 

Rl 
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The existence of the limit ( 11911) — a nontrivial question but we are going to show that 
in all the cases we are interested the limit (I19ip exists. 

We will say that the function / on M is periodical if Vx, y G M such that 7r(x) = 7r(y) 
f(x) = f{y). All the lifts of the functions given on M on M are periodical and all the 
periodical functions can be obtained this way. All the interesting for us dynamic variables 
on M appear exactly as lifts of the functions given on M and that is why are periodical. 
Further the generalized microcanonical distributions on M, built by H(\), Ki(X), ...,Kk(X) 
are also periodical. We are going to demonstrate that if f(x) and w(x) are periodical 
functions then the limit (11 9 1| exists. 

Let us assume that dT 1 is the element of the phase volume on M\. Let us assume that 
Mk[,...,k' is the surface of the level on M, corresponding to the values K[,...,K' k of the 
integrals Kx,...,K^. Let us assume that dv = dT^dipx, dfk is the measure on Mk[,...,k' ■ 
The projections reflects Mr[ k' on ^k' 1 ,...,k' is the surface of the level on M corresponding 
to the values K[, K' k of the integrals K\, Ky.. On Mjc[,...,k' we have a measure dv, That 
can be characterized as a quotient of the phase volume dT by dK\...dKk'- dT = dvdK\...dK^. 
It is evident that the measure dv is a lift of the measure dv from Mk , 1 ,...,k i to Mk[,...,k' by 
virtue of the canonical projection n. To prove the conclusion about the existence of the limit 
(11911) . it is evidently enough to prove that for any K[, K' k the following limit 

lim-1 f f(x)du(x) (192) 

k 



L— >oo 



exists for any (good enough) periodical function f(x) — f o n(x). And for this it is enough, 
evidently, to prove that for any different iable function if) on Mi with a compact support 
there is a limit: 

L^ojk J ^ x ) dT l J ffaiPU-iVk)' (193) 
Mi \<Pi\<L 

But on second thought, the existence of the last limit succeeds directly from k- dimensional 
analogue of the von Neumann ergodic theorem [14], applied to the commutated flows on 
Mk' 1i ...,k' j that are induced by the Hamiltonians K\, ...,K^ keeping the measure dv and to 
restriction of the function / on Mk{,...,k' ■ 

Now in order to establish all the necessary to us thermodynamic relations, for example 
the relation dE = TdS — PdV , we derive them for M cut of Rl, and then tend L to 
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infinity. The only difficulty that can occur is the divergence of the entropy at L — > oo, but 
this divergence is removed by the diminution of the logarithmically divergent constant^ In L 
from the entropy and this constant is the same at any value of the parameter A. 

The inequality (11241) that we need is now derived as before but at A = 0. It is enough to 
obtain again all our conclusions concerning the Ling cell. 

8 Appendix 2. Discussion of the boundary condition 
[T271 

The aim of this appendix is to make more or less obvious the existence of such temperature 
T , above which (SS)e = const. 

Let us assume that at s > s + 5s the integrals becomes inactive. The entropy of 
the system S(E, K') at fixed energy E depends on the value K' of this integral and the 
values of the other integrals that we will not name here. When discussing the equivalency 
principle we said that S(E, K') should not depend on K', but it would be more precise to say 
that S(E,K') reaches the maximum on some interval [a,f3] of the values K, and decreases 
drastically outside it. It is clear that (SS)e = hi(/3 — a). But this kind of dependence 
S(E, K') from K' is not analytical and is proper to the systems with an infinite number of 
degrees of freedom, i.e. in the thermodynamic limit. 

As an example of non-analyticity within the thermodynamic limit we can see the phase 
transition liquid -H- vapor. Let us see the dependence of the system pressure from its volume 
at temperature lower that critical. When the volume is decreasing from infinity the system 
pressure first begins to increase and then on some interval of volume change will remain 
constant. This interval corresponds to the co-existence of the liquid and vapor phases. 
At further volume decrease the pressure will increase again. So this way this isotherm is 
not describe by the analytical function. When temperature is increasing the analyticity is 
restored. The temperature at which the analyticity is restored is called critical. Above this 
temperature the difference between vapor and liquid disappears. 

Another example of the non-analytical behavior of the thermodynamic values which is 
more close to the Ling cell is borrowed from the van der Waals theory of the second-order 
phase transitions [13]. The second-order phase transition is usually related to some symmetry 
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breaking and the breaking of this symmetry is described by a certain parameter of the order 
77. We will analyze here the simplest case of the real-valued rj and Z 2 -symmetry that only 
changes the sign 77; 77 1— > —77. 

In the neighborhood of the point of phase transition the Taylor approximation of the 
thermodynamic potential $(77, T) accurate to the inappreciable terms with higher orders, 
takes the form: 



where T c — critical temperature and Co, c\ — certain positive constants. 

At T > T c $(77, T) has by T one minimum in zero. At T < T c $(77, T) has one local 
maximum by 77 in zero and two minimums in the points 



We assume that the parameter of the order 77 is additive and that is why is lower that critical 
temperature, the analyzed system is separated into domains and in each of them 77 takes one 
of the two opposite values, corresponding to the minimum of the thermodynamic potential 
of the domain. The resulting value 77 will be bounded between two values 771 and 772. From 
the other side, as the domains are big enough we can neglect the interaction between the 
different domains when calculating the thermodynamic potential and the thermodynamic 
potential will be equal to the minimal value of $(77, T) by 77. 

So this way we come to conclusion that instead of $(77, T) we should use its upper 
envelope $(77, T) (by 77) and at T > T c $(77, T) will have the plateau, i.e. $(77, T) will not be 
the analytic function anymore. 

The recovery of the analyticity at the temperature increase is explained in the following 
way. At the temperatures high enough the kinetic energy of every particle becomes much 
more than the potential energy of the particle interaction and the particles can be considered 
non-interacting. I.e. at high temperatures all the particles that compose the system can be 
considered free and for such system all the characteristics must be analytical. 

The same way we will assume that at the temperature high enough S(E ) K') will be the 
analytical function of the parameter K' . Then it is either constant or reaches the maximum 
at an isolated point. But (SS)e = hi(/3 — a), as it was shown above, decreases with the 
temperature increase, i.e. at temperature high enough S(E, K') can not be constant on the 



$(77,T) = c 77 4 + Cl (T-T c )77 2 , 



(194) 




(195) 
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whole interval of the change of K' . I.e. S(E, K') reaches the maximum at an isolated point 
and there is no more difference between using the generalized microcanonical distribution or 
just the microcanonical distribution. 
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